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Abstract 

Consider  a  plane  wave  incident  on  a  conducting  strip  of  infinite  lengtk  coated  on  one  side 
with  a  thin,  uniform  coating  of  lossless  dielectric,  and  let  the  iucident  magnetic  field  vector  be 
parallel  to  the  edges  of  the  strip.  From  first  principles  we  develop  the  integral  equations  relating 
the  total,  incident,  and  scattered  fields  in  the  dielectric  We  expand  the  fields  in  the  dielectric  slab 
as  the  sum  of  three  plane  waves  (the  ‘‘forced*  wave  and  two  oppositely-directed  surface  wa\es) 
The  phase  constants  in  the  x  and  y  directions  are  known  and  the  wave  amplitudes  are  found  using 
least-squared-error  point-match i ng. 

The  far-zone  scattered  field  is  found  by  letting  the  polarization  and  conduction  currents 
radiate  in  free  space. 

Data  are  presented  showing  the  bistatic  scattering  characteristics  of  the  dielectric-coated 
conducting  strip  as  a  function  of  dielectric  permittivity,  slab  width  and  thickne&j,  and  angle  of 
incidence  Results  show  reasonable  agreement  with  conventional  moment  method  solution  using 
pulse  bases  and  Galerkin  testing  functions. 


TE  SCATTERING  FROM  A  DIELECTRIC  COATED  CONDUCTING  STRIP 

I.  Introduction 

J.l  Background 

Today’#  military  environment  increasingly  requires  tactical  and  strategic  aircraft  to  possess 
some  measure  of  “stealth,”  or  low-observable,  capability  The  logic  is  intuitive  if  an  aircraft  cannot 
be  detected,  then  it  cannot  be  shot  down,  and  hence  is  more  effective  than  a  conventional  aircraft 

The  most  potent  threat  an  aircraft  faces  is  acquisition/tracking  radar  in  concert  with  either 
an  infrared-  or  radar-guided  missile  Accordingly,  most  of  the  emphasis  in  the  low-observables 
community  has  been  aimed  at  reducing  the  radar  cross  section  (RCS)  of  aircraft.  The  RCS  of  an 
aircraft  is  a  measure  of  how  much  electromagnetic  (EM)  energy  is  reflected  (scattered)  back  at 
an  illuminating  radar.  Reducing  an  aircraft’s  RCS  reduces  the  range  at  which  a  radar  can  detect 
it.  Although  actual  RCS  measurements  for  current  military  aircraft  are  classified,  a  reasonable 
estimate  is  that  the  RCS  of  newly-developed  low-observable  aircraft  is  comparable  to  that  of  a 
metal  sphere  that  is  10  centimeters  in  diameter  (very  small  indeed’)  (13  22-23). 

Reducing  an  aircraft’s  RCS  is  best  approached  using  a  “divide  and  conquer”  approach  An 
aircraft  may  be  thought  of  as  a  collection  of  sections,  each  of  which  has  an  RCS  associated  with 
it  The  RCS  of  the  total  aircraft  is  then  the  coherent  combination  of  the  individual  RCS’s  This 
partitioning  has  two  advantages. 

First,  if  one  section  of  an  aircraft  has  an  RCS  that  dominates  the  RCS  of  all  other  sections, 
then  the  most  benefit  can  be  attained  by  concentrating  efforts  to  reducing  the  RCS  of  that  one 
dominant  section. 
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Second,  transforming  the  problem  from  the  analysis  of  an  extremely  complex  structure  such 
as  an  aircraft  to  that  of  a  relatively  simple  shape  greatly  reduces  the  difficulty  of  the  problem  In 
most  cases,  only  a  few  canonical  shapes  (e.g ,  cylinders,  plates,  and  spheroids)  are  necessary  to 
adequately  model  an  aircraft  of  arbitrary  configuration 

In  the  same  manner  that  an  aircraft  is  decomposed  into  canonical  shapes,  the  RCS  of  the 
canonical  shapes  may  be  decomposed  into  various  physical  phenomena,  called  scattering  mecha¬ 
nisms.  In  general,  the  scattering  mechanisms  of  any  object  are  dependent  on  the  size,  shape,  and 
material  properties  of  the  object. 

1.2  Statement  of  the  Problem 

One  canonical  shape  of  considerable  interest  »s  a  conducting  strip  covered  on  one  side  with  a 
thin  material  coating  of  known  electrical  properties  For  example,  such  a  geometry  can  be  used  to 
model  an  aircraft  section  covered  by  radar-absorbing  material  (RAM),  or  even  an  airfoil  covered 
with  ice  (6:3-4) 

The  specific  problem  (14  194-202)  to  be  analyzed  in  this  thesis  is  the  electromagnetic  scat¬ 
tering  from  a  dielectric-coated  conducting  strip  of  infinite  length,  as  shown  in  Figure  1  The  strip 


y 


Figure  1.  Two-Dimensional  Scattering  Geometry 

is  arbitrarily  thin,  has  infinite  conductivity  dr,  and  finite  width  u\  The  dielectric  coat.ng  has  width 
tu,  thickness  h,  permittivity  tg  =  crc o  and  permeability  fig  =  UrPo  The  incident  field  is  a  linearly 
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polarized,  uniform  plane  wave  of  unit  amplitude  having  the  electric  field  vector  transverse  to  the 
z-axis,  i.e.,  in  the  transverse  electric  (TE)  sense  (E  in  the  xy-plane,  H  aligned  with  the  z-axis) 

l.S  Approach 

The  first  task  consists  of  reviewing  the  available  literature  concerning  electromagnetic  scat- 
tering  from  coated  objects  The  results  are  compiled  in  Chapter  II. 

The  second  task  consists  of  developing  the  integral  equations  relating  the  scattered  field  to 
the  incident  field  The  total  field  quantities  must  satisfy  the  Helmholtz  equation  as  well  as  the 
classical  boundary  conditions  of  Ej  =  0  at  y  =  and  Ej  rind  Hf  continuous  across  y  =  £  (Fig 
1)  The  results  comprise  Chapter  III. 

The  third  task  is  to  develop  a  moment  method  solution  to  the  integral  equations  using  physical 
bas'13  functions.  The  physical  basis  functions  will  consist  of  a  forced  wave  and  forward  and  reverse 
x-directed  surface  wave  components  (19).  Testing  functions  will  be  oversampled  point-matching 
Dirac  delta  functions,  creating  an  over-determined  linear  system  of  equations  and  a  non-square 
impedance  matrix.  Computer  algorithms  will  be  developed  to  generate  the  impedance  matrix 
elements,  which  will  be  line  and  surface  dimensional  integrals  involving  Hankel  functions  of  the 
second  kind  of  order  0  The  integrals  will  contain  singularities  which  must  be  evaluated  in  the 
principal-value  sense  The  over-determined  system  will  be  solved  in  the  least-square-error  tense 
using  matrix  pseudo-inversion  techniques.  The  results  may  be  found  in  Chapter  IV. 

The  fourth  task  is  to  compare  the  output  of  the  moment  method  algorithms  developed  in 
the  second  task  to  the  output  from  a  standard  moment  method  program  using  pulse  basis/point¬ 
matching.  The  criteria  for  comparison  will  be  absolute  accuracy  of  the  solution,  computer  run  time, 
and  computer  storage  requirements  If  physical  basis  functions  work  as  anticipated,  the  accuracies 
will  be  nearly  identical,  but  the  physical  basis  function  solution  will  require  much  less  computer 
time  and  storage  than  the  standard  moment  method  solution. 
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II.  Literature  Review 


Much  work  has  been  done  in  the  areas  of  “high-frequency”  scattering  theory  -  Geometrical 
Optics  (GO),  the  Geometrical  and  Uniform  Theories  of  Diffraction  (GTD  ic  UTD),  Physical  Optics 
(PO),  and  the  Physical  Theory  of  Diffraction  (PTD)  are  the  most  prevalent.  These  methods 
are  called  “high-frequency”  because  they  are  derived  from  asymptotic  solutions  to  the  scattering 
equations  in  the  hr.  t  as  the  wavelength  of  incident  radiation  becomes  small  compared  to  the 
dimensions  of  scaiterer.  In  this  thesis,  however,  the  scatter  is  on  the  order  of  a  wavelength  (a 
few  to  tens  of  wavelengths  in  width)  so  the  “high-frequency”  methods  are  uot  strictly  applicable 
Therefore,  only  material  dealing  with  the  method  of  moments  implementation  will  be  reviewed 
The  review  divides  into  two  distinct  portions  material  dealing  with  developing  the  exact  scattering 
integral  equations,  and  that  dealing  with  approximating  the  solution  to  the  integral  equation. 

2.1  Integral  Equations 

The  integral  equations  relating  the  scattered  e**ctric  field  E'  to  the  incident  electric  field  E* 
may  be  derived  using  several  techniques;  all  are  based  on  fundamental  EM  theory  and  are  exact. 
In  general,  developing  the  integral  equation  is  a  two-step  process 

2.1.1  Equivalent  Currents  As  a  first  step,  we  may  remove  the  conducting  portion  of  the 
scattering  body  and  replace  it  with  fret  space  and  equivalent  electric  and  magnetic  surface  currents, 
J  and  M,  on  the  external  boundaries  of  the  body  (2, 10, 14, 21).  This  is  a  direct  application  of  the 
surface  equivalence  theorem  (3:320-333). 

The  surface  currents  are  such  that  they  radiate  E*  outside  the  conducting  scattering  body 
and  —E*  within  it.  Thu  creates  (necessarily)  a  vanishing  field  inside  the  conducting  body,  the 
technique  is  sometimes  referred  to  as  the  “null  field”  approach  (28)  In  addition,  the  fields  radiated 
by  the  surface  currents  must  satiny  the  appropriate  boundary  conditions: 
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1.  Continuity  of  the  tangential  electric  and  magnetic  fields  at  a  dielectric  interface. 

2.  Vanishing  tangential  electric  field  at  the  surface  of  a  perfect  conductor. 

In  a  similar  manner,  we  may  replace  the  dielectric  portion  of  the  scatterer  by  free  space  and 
equivalent  volumetric  current  densities.  These  volumetric  currents  must  satisfy  the  same  boundary 
conditions  as  the  surface  currents.  Once  the  surface  and  volume  currents  are  specified,  the  scattered 
fields  £'  may  be  easily  found  (21,  22).  The  problem  is  thus  converted  into  finding  the  equivalent 
surface  currents  J  and  M  in  terms  of  the  incident  field  E‘ 

2  1.2  Integral  Equation  Formulation  The  second  step  is  to  formulate  the  integral  equations 
for  equivalent  currents  J  and  Af .  There  are  many  ways  to  write  the  integral  equations;  the  method 
used  is  influenced  by  the  nature  of  the  scattering  body,  the  polarization  of  the  incident  field  E' ,  or 
even  the  personal  whims  of  the  author. 

Kent  (14)  uses  Rumsey’t  Reaction  Integral  Equation  for  the  specific  cases  of  a  finite-length  di¬ 
electric-coated  wire,  an  infinite  dielectric  slab,  and  a  dielectric-coated  conducting  strip.  Rogers  (21) 
uses  the  dyadic  free-space  Green’s  function  to  write  an  integral  equation  defining  the  surface  cur¬ 
rents  J  and  M.  Use  of  dyadic®  is  dictated  by  the  arbitrary  polarization  of  the  incident  field  E'. 
Rogers  asserts  that  the  resulting  equation  is  general  and  can  be  manipulated  into  Mitschang’s  vector 
and  scalar  potential  function,  the  electric  field  integral  equation  (EFIE),  the  magnetic  field  integral 
equation  (MFIE),  and  ihe  combined  field  integral  equation  (CFIE).  Huddleston  (10)  uses  a  CF1E 
formulation  in  which  the  scattered  field  E1  is  expressed  in  terms  of  integrodifferential  operators 
involving  dyadic  free-space  Green’s  functions.  Roj?s  (22)  examines  the  inhomogeneous  dielectric 
cylindrical  scatterer;  he  uses  the  volume  equivalence  theorem  (3:327-329)  to  define  volumetric  cur¬ 
rent  densities  within  the  cyT.ader  in  a  manner  analogous  to  that  involving  surface  current  densities. 
Again,  Rojas’  solution  uses  the  dyadic  free-space  Green’s  function  in  the  final  integral  equation. 
Arvas  (2)  uses  the  standard  EFIE  solution  for  the  scattering  from  multiple  perfectly  conducting  and 
lousy  dielectric  cylinders  of  arbitrary  cross-section,  including  dielectric-coated  conducting  cylinders 


2.2  Method  of  Moments 


Once  the  integral  equation  is  written,  we  most  solve  it  to  find  the  scattered  field  E*  in  terms 
of  the  incident  field  E*  .  As  stated  earlier,  no  analytic,  closed-form  solution  exists  in  general, 
and  we  mus*  'inloy  numerical  techniques  to  approximate  a  solution.  The  preferred  technique, 
called  the  Method  of  Moments  (9),  is  to  transform  the  integral  equation  into  a  system  of  linear 
algebraic  equations  This  is  done  by  expanding  the  unknown  currents  in  basis  functions  and  then 
introducing  testing  functions.  Commonly,  the  same  expansion  functions  are  used  as  both  basis  and 
testing  functions;  this  is  known  as  Gaierkin’s  method. 

2.2.1  Subsecitonal  Basis  Functions  Basis  functions  which  exist  only  on  discrete  positional 
subsections  are  known  as  subsectional  basis  functions  (3-683)  Examples  include  pulses  and  piece- 
wise  sinusoids  for  two-dimensional  geometries,  and  square  and  triangular  patches  for  three-dimen¬ 
sional  geometries.  Subsectional  basis  functions  provide  the  flexibility  necessary  to  analyze  arbitrary 
geometries  where  the  nature  of  J  and  M  is  not  known.  However,  solutions  using  subsectional  basis 
functions  typically  require  a  large  number  of  them  to  achieve  a  given  level  of  accuracy  in  the 
solution. 

Huddleston  (10)  uses  two-dimensional  triangular  basis  functions  and  GrJerkin  testing  func¬ 
tions  to  expand  the  surface  currents  for  coated  and  uncoated  spheres,  closed  circular  cylinders,  and 
flat-back  circular  cones.  Rogers  (21)  uses  bi-triangular  basis  functions  for  the  case  of  coated  and 
uncoated  fiat  square  plates  Arvas  (2)  and  Rojas  (22)  use  pulse  basis/point  matching  to  analyze 
several  varied  two-dimensional  geometries  involving  dielectric  and/or  conducting  cylinders 

Most  authors  are  eager  to  report  agreement  with  measured  and/or  other  published  data  but 
seem  less  enthusiastic  to  discuss  the  number  of  sub-domain  basis  functions  N  or  the  computer  time 
and  storage  requirements  necessary  to  achieve  their  results.  Only  Rojas  (22)  states  the  number  of 
expansion  functions  N  and  the  computer  run-time  needed  to  run  his  code.  For  example,  he  relates 
that  for  a  circular-shell  dielectric  cylinder,  the  system  of  linear  equations  to  be  solved  is  of  order 


N  =  246  and  the  computer  code  requires  3  minutes  to  run  on  a  VAX  11/780. 


2.2.2  Enitrc-Domatn  Basis  Functions  Basis  functions  which  expand  the  unknown  surface 
currents  J  and  M  into  “modes”  which  exist  over  the  entire  scattering  body  are  known  as  entire- 
domain  basis  functions  Examples  include  Fourier  aeries  and  Chebysfcev  polynomials  Entire- 
domain  basis  functions  require  knowledge  of  the  nature  of  unknown  J  and  Af,  but  can  greatly 
reduce  the  order  N  of  the  system  to  be  solved  (3.685). 

Kent  (14)  uses  a  special  form  of  the  entire-domain  basis  function  called  the  physical  basis 
function  (PBF),  introduced  by  Richmond  in  1983  (19).  PBFs  can  be  used  to 

analyze  particular  structures  which  may  be  large  in  one  dimension  but  which  will  require 
only  a  small  number  of  expansion  functions  for  the  unknown  currents.  The  method  will 
have  the  additional  property  that  most  or  all  of  the  current  components  have  specific 
physical  interpretations. . . 

Specifically,  Kent  models  the  equivalent  polarization  and  conduction  currents 

as  being  composed  of  a  ‘forced’  wave  (which  is  similar  to  the  classical  physical  optics 
(PO)  current)  and  two  oppositely  directed  surface  waves. 

Using  PBFs  and  Galerkin  testing  functions,  the  order  of  the  system  of  equations  is  N=3,  a  very 
small  number  indeed  compared  to  subsectional  basis  functions.  The  price  paid  is  the  increased 
complexity  of  the  numerical  integrations  which  must  be  performed  to  fill  the  3x3  impedance 
matrix. 

For  the  dielectric-coated  wire,  Kent  states  that  the  PBF  solution  agrees  well  with  the  standard 
moment  method  solution  using  subsectional  piecewise-sinusoidal  (PWS)  expansion  functions  for 
scatterer  length  greater  than  two  wavelengths  Some  disagreement  exists  for  shorter  lengths;  Kent 
suggests  that  a  traveling  wave  current  exists  on  the  structure  which  is  not  accounted  for  in  the 
PBF  solution. 
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For  the  dielectric  slab,  the  author  report®  that  the  PBF  solution  shows  excellent  agreement 
with  the  PWS  solution.  Kent  notes  that  the  reverse  surface  wave  component  provides  the  dominant 
scattered  field  E*  at  near-grazing  angles 

2.2  3  Hybrid  Basis  Functions  Several  attempts  have  been  made  to  use  basis  functions  that 
take  the  form  of  asymptotic  solutions  to  high-frequency  scattering  problems.  The  goal  is  to  reduce 
the  number  of  unknowns  sufficiently  so  that  the  method  of  moments  can  be  applied  to  electrically 
large  scatterers.  This  technique  seems  limited  to  perfectly  electrically  conducting  (PEC)  scatterers 

Blue  (5)  examines  a  two-dimensional  (PEC)  strip.  He  hypothesizes  that  the  strip  currents 
may  be  well  approximated  by  physical  optics  in  the  center  of  the  strip,  and  by  the  Geometrical 
Theory  of  Diffraction  (GTD)  near  the  edges.  He  thus  expands  the  unknown  currents  as  several 
physical-optics  components  and  several  oppositely-directed  diffraction  components  He  allows  the 
PO  currents  to  exist  only  away  from  the  edges,  but  allows  the  GTD  currents  to  exist  everywhere 
on  the  strip.  He  then  solves  for  the  unknown  amplitudes  of  these  components  by  forcing  the  total 
current  to  satisfy  (in  a  least-squared-error  sense)  the  integral  equation  at  M  match  points,  where 
Af  is  greater  than  or  equal  to  the  number  of  unknown  coefficients.  For  a  PEC  strip  that  is  200 
wavelengths  wide,  Blue  achieves  good  results  with  just  17  basis  functions  and  17  match  points. 

Another  hybrid  technique  involves  using  GTD  currents  where  the  scatterer  is  smooth  and 
subdomain  pulses  near  edges  and  other  discontinuities.  Burnside  (7)  analyzes  a  2A  PEC  square 
cylinder  using  16  pulses  at  the  four  corners  and  eight  diffracted  currents,  for  a  total  of  24  unknowns 
The  unknowns  are  found  using  least-squares  point  matching.  The  results  agree  very  well  with  a 
standard  moment  methods  using  straight  pulse/point-matching  with  56  unknowns. 

Thiele  and  Newhouse  (26)  modify  the  generalized  impedance  matrix  to  include  contributions 
from  GTD  currents,  with  an  emphasis  on  near-field  interactions.  Specifically,  they  investigate  a 
monopole  on  a  PEC  wedge,  a  monopole  on  a  PEC  circular  disc,  and  a  monopole  near  a  PEC  step 
Inclusion  of  GTD  effects  allows  Thiele  and  Newhaouae  to  reduce  the  size  of  the  impedance  matrix 
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while  maintaining  good  agreement  with  measured  data  They  suggests  that 

..  a  potential  advantage  of  such  an  approach  is  that  in  many  problems  the  impedance 

matrix  may  be  small  enough  to  calculate  data  at  sufficiently  many  frequencies  to  obtain 

time  domain  solutions  via  Fourier  transformation... 

Sahalos  and  Thiele  (23)  extend  the  techniques  to  three  dimensions  by  analyzing  the  scattering 
from  a  three-sided  PEC  pyramid  After  achieving  encouraging  results,  the  authors  suggest  that  the 
GTD-MM  technique  may  be  applied  to  a  PEC  pyramid  with  a  dielectric  coating  on  one  or  more 
faces. 

The  hybrid  techniques  are  similar  to  physical  basis  functions  in  that  both  approaches  attempt 
to  exploit  the  expected  form  of  the  scattered  fields.  The  two  approaches  differ  in  two  important 
respects,  however.  First,  hybrid  techniques  use  functions  from  the  asymptotic  high-frequency  so¬ 
lutions  of  GTD  and  PO  while  physical  basis  functions  use  functions  from  the  corresponding  one¬ 
dimensional  scatterer.  Second,  hybrid  techniques  tend  to  be  subsectional  m  that  they  use  different 
basis  functions  over  different  parts  of  the  scatterer  based  on  proximity  to  discontinuities.  On  the 
other  hand,  physical  basis  functions  are  entire  basis  functions  and  exist  over  the  entire  scatterer, 
even  near  the  edges. 

£.3  Implications  to  Problem  at  Hand 

Much  work  has  been  done  in  the  analysis  of  the  low-frequency  EM  scattering  from  two- 
dimensional  dielectric-coated  conducting  cylinders  of  arbitrary  cross-section.  The  standard  method 
is  to  replace  the  unknown  scattered  field  E*  with  unknown  equivalent  surface  and  volumetric  current 
densities  J  and  M ,  write  integral  equations  for  J  and  Af,  and  find  approximate  solutions  to  the 
integral  equations  using  the  moment  method. 

The  dielectric-coated  conducting  strip  is  a  specific  case  of  the  dielectric-coated  conducting 
cylinder  of  arbitrary  cross-section.  It  may  be  replaced  by  the  equivalent  surface  and  volume  current 
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problem,  and  integral  equations  written  to  relate  J  and  M  to  the  incident  field  E*.  We  suspect 
that  currents  are  dominated  by  the  “forced”  wave  and  two  surface  waves  that  comprise  the  physical 
basis  functions,  and  that  a  moment  method  solution  using  PBFs  should  yield  an  accurate  solution 
while  minimizing  computer  time  and  storage  requirements 
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III.  Theory 


3  J  Scattering  Geometry 

Consider  a  perfectly  conducting  thin  strip  of  finite  width  w  and  infinite  length,  coated  on 
one  side  with  a  material  slab  of  width  u>  and  thickness  h ,  as  shown  in  Figure  1.  The  material 
is  characterized  by  permeability  fid  =  firfio  and  permittivity  a  =  erco,  where  fio  and  Co  are 
the  permeability  and  permittivity  of  free  space,  respectively,  fir  and  cr  are,  in  general,  complex 
quantities  with  real  parts  greater  than  unity.  If  fir  and  er  are  real,  then  the  material  is  lossless 

The  time-harmonic  incident  wave  is  a  transverse-electric  (TE)  polarized  uniform  plane  wave, 
with  the  electric  field  vector  E  in  the  xy-plane  and  the  magnetic  field  vector  H  parallel  to  the 
z-axis.  The  time  dependence  efut  is  assumed  and  suppressed. 

3.2  Integral  Equaitons 

The  total  electric  and  magnetic  fields  are  composed  of  incident  and  scattered  components 


r 

=  E‘  +  E* 

t  V/m] 

(1) 

Hr 

=  H'  +  H* 

[A/m] 

(2) 

We  define  the  incident  field  (E*  and  H*)  to  be  that  which  would  exist  if  the  scaiterer  were  not 
present  and  the  total  field  (Er  and  Hr)  to  be  that  which  exists  in  the  presence  of  the  scatterer. 
The  scattered  field  (E*  and  H*)  is  the  thus  the  difference  between  the  incident  and  total  fields. 

All  fields  must  satisfy  Maxwell’s  equations  for  time-harmonic  electromagnetic  waves.  Addi¬ 
tionally,  the  scattering  geometry  dictates  that  certain  boundary  conditions  must  be  met.  The  first 
boundary  condition  requires  the  total  tangential  electric  field  to  vanish  at  the  conductor  (8:34),  or 
mathematically,  n  x  Er  =  0  for  \x\  <  *  and  y  =  —  where  n  is  the  unit  vector  normal  to  the 
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conductor  surface.  The  second  boundary  condition  requires  continuity  of  the  tangential  electric 
and  magnetic  fields  across  the  material/frec-space  interface  (3*14-15) 

S.2.1  Equivalent  Currents  The  surface  and  volume  equivalence  theorems  (3,  8)  allow  us  to 
replace  the  scatterer  with  free  space  and  equivalent  electric  and  magnetic  currents  which  radiate 
the  scattered  fields  as  shown  in  Figure  2.  A  conduction  surface  cunent  density  (8:34) 


Jc  =  S  Jl 

=  u  x  H3"  =  y  x  i.H'f  =  xHj  [d/m] 


(3) 


exists  on  the  conductor  (|r|  <  =  -3)  while  equivalent  polarization  electric  and  magnetic 

current  densities  (3*328) 

=  9J?  +  yJ? 

=  juc.(cr -l)ET  [d/m2) 

m*<  =  iA/;» 

=  J<j(i.(p,  -  1)Ht  [v/m’j 


(4) 

(5) 
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exist  within  the  material  (|rj  <  5*  Ivl  <  f)  Notice  the  difference  in  units  between  the  conduction 
and  equivalent  current  densities.  Je ,  J<f,  and  M<f  radiate  the  scattered  fields  E*  and  H'  in  free 
space. 

3.2.2  Vector  Potentials  We  may  define  a  magnetic  vector  potential  A  as  the  integration  of 
the  electric  current  density  J  against  a  kernel  function  (8:80): 


-  f  1 1  - 

A«  =  T&ZFi***  <6> 

r  =  xx  +  yy  +  az  (7) 

|r-F|  =  \J(x  -  z’Y  +  (y  -  y)3  +  (z  -  z'f  (8) 

to  =  u,/»<oto  (9) 


The  scatterer  is  two-dimensional,  so  J  is  not  a  function  of  s'  (see  Eqs  3  and  4)  and  the 
integration  in  s'  may  be  evaluated  analytically  using  the  identity  (3:581) 

(10> 

where  H ^  is  the  Hankel  function  of  the  second  kind  of  order  0  Furthermore,  the  currents  exist 
only  on  the  strip  and  within  the  material  (|x|  <  f  and  |$/|  <  £),  so  the  magnetic  vector  potential 
is  then 

A(x,v)  =  j*f*  S(x'.  J)  (ioV'(r-r')J  +  (!/-y')5)  dz'dj  (11) 

The  integration  kernel 

lf/<s)(W(*-*')s  +  (v-  7y)  (i2) 

is  the  two-dimensional  free-space  Green’s  function  (3.902),  G(z,  jr,x',j/),  so  we  may  equivalently 
write 
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A(z,y)  =  f  f  3{z',  ]/)G(zty,z',t/)dx'dj/ 
J-V-* 


The  electric  current  density  J  is  composed  of  two  parts,  the  one-dimensional  conduction 
current  density  3C  (Equation  3)  and  the  two-dimensional  equivalent  polarization  current  density 
Je#  (Equation  4).  Thus,  by  superposition,  the  magnetic  vector  potential  may  be  written  as  the 
sum  of  two  integrals 

A(x,v)  =  J  ^3%z‘)G(z,y,z',- |)*t' 

+  J*J^3’'(z',tf)G{z,r,z',S)dz'<W  (14) 

Notice  that  j/  =  —  |  in  the  line  integral  of  Equation  14. 

By  duality  (8:98-100),  an  electric  vector  potential  F  can  be  defined  as  a  function  of  the 
equivalent  magnetic  current  density  M  =  Mff 

.  ft  f *  . . . . 


F(r,y)=  /  /  M,f  (*', }/) G(zt yt x\ tf) dz* d]/ 


Notice  that  there  is  no  one-dimensional  magnetic  conduction  current  on  the  strip. 


3  S  3  Scattered  Fields  in  Terms  of  Vector  Potentials  The  scattered  fields  (E*  and  H*)  are 
related  to  the  vector  potentials  (A*  and  F*)  and  current  densities  (J  and  M)  by  (8:99) 


-V  x  F  +  - — (y  x  V  x  A  —  J) 
JUCq 

VxA+  (V  xyxF-M) 
JUflQ 


Equations  16  and  17  may  be  simplified  by  noting  that  the  electric  and  magnetic  vector  potentials 
satisfy  the  inhomogeneous  Helmholtz  equations  (8:77) 


Vs  A  +  tjA  =  -J 
y’F  +  tjF  =  -M 
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and  making  use  of  the  vector  identity  ^xvxi?=V(V  >?)  -  vV  (3-927)  We  may  then  recast 
Equations  16  and  17  as 


E'  =  -VxE-yw/ioA  +  T^—  V(V  •  A) 

JUCo 

H*  =  V  x  A  -  jue0F  +  -A- V(V  F) 
JU)to 


(20) 

(21) 


To  make  use  of  Equations  20  and  21,  we  will  need  to  define  the  the  gradient  of  the  divergence 
(V[V  •  </*))  and  the  curl  (V  x  tO  operators  acting  on  A  and  F  Noting  that  £  =  0,  the  elementary 
operators  (gradient,  divergence,  and  curl)  acting  on  arbitrary  scalar  and  vector  functions  V  and 
c!1  ■=  v'i x  +  Vj-y  +  t M,  respectively,  are  given  by  (3:925) 


V-t? 


a^ 

ax  yay 


ax 

.w. 

8y 


a» 


-  * 


at>, 

ax 


+  t 


X  ax  “  ay/ 


(22) 

(23) 

(24) 


We  construct  the  gradient  of  the  divergence  operator  using  Equations  22  and  23  as 


V(V  ■  W 


U±  +  y±)  (9Jl.  + 

V  8x  y8y/  v  8x  +  ay  ) 

ra^,  a^l 

l  ax’  axayj  +  y 


fay, 

\8x8y 


(25) 


$.£.4  Vector  Potential  in  Term*  of  Cxrrept  Densities  It  b  convenient  to  think  of  the  vector 
potentials  as  being  composed  of  three  components,  each  radiated  by  one  of  the  three  current 
densities  Je,  J'f,  or  M<#.  Then 


A‘(*,y)  =  iA;(x,y) 

=  i[^J',(x')G{x,y,x',~)dx'  (26) 

A'*(x,y)  =  fMr(x.y)  +  JAJ«(x,y) 
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(27) 


=  /*  /*  {*W> l/)  +  G{z%yt*?ti/Wd)/ 

F«(zty)  =  zF?(xty) 

=  *  r  /*  Affix' ttf)G(z,y,zW)dz'dj/  (28) 

In  Equations  26, 27,  and  28,  the  current  densities  are  functions  of  primed  variables,  and  the  Green’s 
function,  G,  (see  Equation  12)  is  a  function  of  pnmed  and  unprimed  coordinates  In  the  following 
development,  the  arguments  of  the  current  densities  and  G  are  suppressed  for  brevity. 

To  use  Equations  20  and  21,  we  need  the  curl  and  the  divergence  of  the  gradient  of  the  three 
vector  potential  components.  Noting  that  =:  A\  ss  A]1  =  =  F^f  =  0  and  using  Equations 

24,  25,  26,  27,  and  28,  we  write 


VXA  =  *{±A<  +  ±AV  -  ±A‘<}  (29) 

*xf ^  (30) 

V(VA)  =  *{|U  +  +  jg^} 

+  mA‘ +  wM  (31) 

V(VF)  =  0  (32) 


We  may  now  substitute  Equations  26, 27, and  28  into  Equations  20-32  to  yield 

(33) 

=  {4y-  *£}  J.\  I.\  <U) 

v(vX)  = 

+  &££***«} 
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f{ili *  A/.'/.',':'0'-'" 


V(V-F)  =  0 


(35) 

(36) 


3.2. 5  The  Scattered  Fields  Equations  33-36  can  be  decomposed  into  parts  associated  with 
J',  J*f,  7**,  and  Af/f.  We  can  then  find  the  scattered  electric  field  E*  as  a  function  of  J£,  J*1, 
J*1,  and  Af/f  using  Equations  20,  34,  14,  and  35.  This  gives 

E’  =  E'(Jx)  +  E'(J*f)  +  +  E*(A//«)  (37) 

E'(J')  =  -jvyoxf*  JIGdx' 

+M*&L\J‘Gd*' +  *£ryL\J‘Gdz'}  (38) 

-  /■!  rf 

E'W)  =  J  Ji'Giz'dtf 

E*(j;»)  =  -ju/ioy  j  ^  J*1  Gdx'dt/ 

1  f  85  /I  8s  /I  7?  ' 

E  •(«;«)  =  — Ml'Gdi'd-J  +  yA  j'j^WGdx'd;/  (41) 

(42) 


Equivalent  expressions  for  the  scattered  magnetic  field  H*  are  considerably  simpler  since 
y(V  •  F)  =  0,  and  may  be  found  using  Equations  21,  33, 15,  and  36 


*«*>  =  -ifJJ^i'cdz'w 


(43) 

(44) 
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H '(«/*)  =  -) uc0z  /  /  M'.'Gdz'd ]/ 


(45) 

(46) 


Wc  can  restructure  Equations  38-41  to  find  the  i  and  y  components  of  E*  Noting  that 
kl  =  w2poco»  we  find 


=  jk{ll  +  &LffJ'Gdx' 

E-w  =  +  &}L\L\j',Gd*,d'/ 

W)  =  7~mL\L\j;,GdI'di' 
BUM!')  =  -jjj*  J*  UPG&# 

e;(js)  =  t1-/!,  r  Jicdx' 

9  jutQ  ozdy 

1  ^  r? 

Wf)  =  /  VG&& 

*  *  jweo&tfy 


(47) 

(48) 

(49) 

(50) 

(51) 

(52) 

(53) 

(54) 

(55) 


TAe  Final  Integral  Equations  The  material  coating  of  interest  in  this  thesis  is  a  non¬ 
magnetic  dielectric.  Then  &  =  1  and  by  Equation  5  no  equivalent  magnetic  current  M  exists  and 
the  scattered  fields  ere  given  by  Equations  43-45  (for  H*),  48-50  (for  E*),  and  52-54  (for  £J)  We 
now  wish  to  write  a  set  of  integral  equations  in  which  the  only  unknowns  are  the  current  densities 
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and  Jc.  This  set  is 


p. 


p. 


juco(t, 

+ 


1  1  (  6*  ft 

-  TJro{e^UUJi,Gdt'^ 

+ [kl+W-J-\jr  ciz'd< +  &L\j‘Gdi} 
H-  +  £I!£*g*« 

8  rt  Q  ft 

-*jJ.fVG**f  +  eiU*0* 


(56) 


(57) 


(58) 


Equations  56-58  are  three  coupled  integral  equations  with  unknown  quantities  J**,  J*f ,  and 
J*,  and  agree  with  Kent’s  development  (15).  In  Chapter  IV,  we  apply  the  method  of  moments  to 
find  the  unknown  currents. 


3.3  Two-Dimensional  Scattering  Width 

We  are  interested  iu  finding  the  two-dimensional  scattering  width  C2D  of  the  dielectric-coated 
conductings*  rip,  both  monostatically  (in  the  direction  of  the  incident  wave)  and  bistatically  (not  in 
the  direction  of  the  incident  wave),  as  shown  in  Figure  3  The  scattering  width  is  given  by  (3.577) 


<T2D  “  lim 
^—•00 


ixP-Z-j 

.  |E‘I 


(59) 


If  we  let  the  incident  plane  wave  have  unit  amplitude,  |E*|  =  1,  then  Equation  59  becomes  simply 


<?2D  =  lim  2irpJE*|2  (60) 

^*co 

Using  Equations  47-55,  and  assuming  the  observation  point  (xrp)  hes  outside  the  scatterer, 
we  may  move  the  derivatives  inside  the  integrals  (This  is  valid  since  the  singularity  of  the  Green’s 
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function  is  not  in  the  region  of  integration.)  Substituting  expressions  for  and  (see 

Appendix  S),  we  find  the  scattered  electric  field  is  given  by 


dz' 


+ ["‘mt5 ■ dz'd!/ 


f  j;,(x'.v')(l  ~  *'$  ~  y')/ff)(*„fi)<fcw}  (61) 

$(.,*)  =  ^  { /_*  W~ 

* r/I  w*  ps?  -  **}<"> 


where  R  =  ^/(r  -  x')2  +  (y  -  y')2. 


We  are  interested  in  the  scattered  field  far  from  the  scatterer,  i.e.,  in  the  limit  as  p  -♦  oo. 
We  may  then  replace  the  Hankel  functions  by  their  large- argument  asymptotic  forms  (17.364).  For 
amplitude  terms,  we  may  replace  R  with  p,  (r  —  z')  with  p coetf,  and  (y  -  y*)  with  p sintf.  For 
phase  terms,  we  may  replace  R  with  p  —  r*cos^  -  y'sin^.  We  also  see  that  as  p  —►  oo,  the  terms 
containing  }j\7\koR)  become  subdominant  to  terms  containing  H^(koR)  and  may  be  neglected 
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We  may  then  write 


£;(z,v)  =  ^  /III  gm^/ 

4  V  P 

(63) 

EUz.y)  = 

4  y 

(64) 

+  «in<S  t‘***<~*W"‘*l&  dtf 

-cos(s  yy*  »>)  ejM*,co.*+»,«o«(,I<</^| 

(65) 

If  we  express  the  scattered  electric  field  in  cylindrical  coordinates  (3  922),  we  find 

E J  =  cos  <j>E'g  +  sin^EJ  =0  (66) 

=  —sin  <j>E'g  +  cosjE* 

•  1  <”> 

As  expected,  the  scattered  field  has  no  radial  component  in  the  far  field 

Plugging  Equation  67  into  Equation  60,  we  now  write  the  two-dimensional  scattering  width 
as 

|/p  (68) 

Once  the  currents  JJ,  J*f ,  and  J*f  are  found,  Equations  65  and  68  may  be  used  to  find  the 
scattering  width  of  the  dielectric-coated  conducting  strip  for  any  observation  angle  <£. 
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IV.  Moment  Method  Solution  Using  Physical  Basis  Functions 


Equations  56-58  define  the  unknown  current  densities  J**,  J*1,  and  in  terms  of  the  known 
incident  field  components  ££,  and  H\.  Unfortunately,  these  integral  equations  cannot  be  solved 
analytically.  We  apply  the  method  of  moments  (9)  and  transform  the  integral  equations  into  a  set 
of  linear  algebraic  equations  which  can  be  solved  numerically. 


4J  Basts  Functions 

The  first  step  in  the  method  of  moments  is  to  select  a  set  of  basis  functions  in  which  to  expand 
the  unknown  quantities. 


•/;*(*.  v) 
•/;•(*.  y) 


]£  On /„*(*.») 


E‘n/n‘(x,V) 


(69) 

(70) 

(71) 


We  choose  to  use  “physical  basis  functions”  as  defined  by  Richmond  (19,  20)  and  Kent  (14)  The 
physical  basis  functions  consist  of  a  “forced”  wave  and  two  oppositely-directed  surface  waves  (25) 
The  forced  wave  is  the  field  that  would  exist  in  the  dielectric  if  it  were  infinitely  wide  (u»  -♦  00) 
The  surface  waves  are  the  dominant  zero-order  TM1  surface  waves  which  have  a  cutoff  frequency 
of  zero.  The  forced  wave  and  surface  waves  are  developed  in  detail  in  Appendix  A  In  terms  of  the 
three  physical  basis  functions,  the  field  and  current  components  in  the  dielectric  may  be  written  as 


E‘{ >•*>  =  hcZ'-l)  =  ti  An  C"  ’iD  [/n  (V+  81  ' 


^(x’v) = iSA) = £  lfn  (w+ 9] 


#f(*.V)  =  53  Anc,n  coe  |/„  (t /  +  trSl‘‘ 
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£(*)  =  H?  (*.  ■ -|)  =  £  -4,  c,„  «-»•*  (75) 

The  terms  An  describe  the  amplitudes  of  the  physical  basis  functions  relative  to  each  other  The 
“field  amplitudes”  crn,  cy„»  and  c,n  describe  the  relative  amplitudes  of  the  individual  field  compo- 
nents  of  the  physical  basis  functions.  The  field  amplitudes  are  defined  by  the  form  of  the  physical 
basis  functions;  they  are  developed  in  Appendix  A  and  given  by 

Csl  cri  cx  s 
Cyl  Cy  2  Cy3 
Cil  C,2  C,3 

An  important  property  of  the  physical  basis  functions,  and  one  that  sets  them  apart  from  pulse 
bases,  piecewise-sinusoids,  etc.,  is  that  they  sene  to  couple  the  x-,  y~  and  x-directed  field  com¬ 
ponents  Because  the  physical  basis  functions  individually  satisfy  Maxwell’s  equations,  if  we  know 
then  we  automatically  know  and  Hj.  So  instead  of  needing  to  solve  for  the  3Af  unknowns 
of  Equations  69-  71  (alt  02, . . . ,  ctt),  we  need  only  solve  for  the  three  wave  amplitudes  An. 

Propagation  constants  fn  and  gn  are  chosen  such  that  n  =  1  represents  the  forced  wave,  n  =  2 
represents  the  4-x-directed  surface  wave,  and  n  —  3  represents  the  -x-directed  surface  wave  /„ 
and  gn  are  real  and  functions  of  the  relative  permittivity  cr  and  thickness  h  of  the  dielectric,  and 
the  angle  of  incidence  8  and  frequency  u  of  the  incident  wave.  Specific  equations  for  f„  and  gn  are 
given  in  Appendix  A;  in  this  chapter  they  are  treated  as  known  constants. 

Expanding  the  currents  in  Equations  56-58,  we  find 

£(*.»)  =  X)  ,in  |/n  (»+ 5)] 

-  (‘r  -  1)  JJ^  Xj  A"  **»  *in  [/»  +  |)]  '* *' G  dz>  W 

■(t'_1)eS;//s 't «*[/»(»’ +  5)] GdI'd,/ 


j  cos  8a  -h—  ill  - 

'  *  JVt4  JUt  A 

sin  $d  rf2 

wc  4  ut  4 

•Jtr/Vo  -1  -1 


(76) 
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H',(z,y)  =  ]Tj  coe  [/»  (»+  |)] 

£  jjs  E  ^ c” sin  ['» (y + 1)]  «■’  '•  *’ c  di>  ^ 

~h  j j  £  ^  c»»  c“  [/-  (*  • +  3)]  '• *' G  dz' v 

+Tyjc±^-‘’''’Gdz‘  <79> 

where  G  =  -  x')2  4-  (y  —  i/)2),  region  S  is  the  cross-sectional  area  of  the  dielectric 

slab  (|x'|  <  *,  \t/\  <  £)  and  contour  C  lies  along  the  conducting  strip  (jx'j  <  f,!/  =  — f ) 

4.2  Testing  Functions  and  Match  Potnls 

In  order  to  obtain  a  system  of  linear  equations  to  solve  for  the  three  unknown  wave  amplitudes, 
we  must  define  a  set  of  testing  functions.  We  choose  to  use  M  (where  M  >3)  Dirac  delta  functions 
which  "fire”  at  the  points  {(x,y)m}^f2l  as  the  testing  functions;  this  is  sometimes  referred  to 
as  “point  matching"  where  the  points  (x,y)m  are  called  “match  points"  This  is  equivalent  to 
enforcing,  in  a  least-squares  sense,  ET(x,y)m  ~  E‘(x,y)m  +  E*(x,y)m,  m  =  1,  M . 

The  M  match  points  may  lie  either  on  the  conductor  or  within  the  dielectric  slab.  No 
guidelines  exist  to  tell  us  where  to  place  the  match  points,  so  we  must  write  the  moment  method 
algorithms  to  allow  many  different  configurations  of  match  points  and  then  look  for  convergence 
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to  a  common  solution.  However,  we  must  enforce  some  order  in  the  match  point  location,  so  we 
partition  the  dielectric  slab  into  a  grid  of  rectangles,  as  shown  in  Figure  4,  with  each  rectangle  of 
size  Ax  x  Ay.  We  pick  Ax,  Ay  «  A j  =  y/s^Xo  such  that  an  integral  number  of  rectangular  cells 


will  fit  inside  the  dielectric.  Then  we  may  write 


1— At— 1 

r _ 

r 

Ay 

JL 

■ 

■M 

yj 

yi 

Z 

1  * 

2 

X 

/ 

Figure  4.  Match  Point  Locations  Withir  Dielectric  Slab 


(/  +  l)Ax 

(80) 

(-M-l)Ay 

(81) 

We  then  pick  IxJ  match  points  m  the  interior  r  f  the  dielectric  and  /  match  points  on  the  conductor, 
for  a  total  of  M  =  /( J  + 1)  match  points.  Wr  may  then  codify  the  location  of  the  mth  match  point 


(*.y)m 

=  (*.,%) 

(82) 

where 

w 

X|  =  --  +  .Ax, 

«=  1 ,/ 

(83) 

h  ,A 

V,  =  ~2+jAy, 

V,. 

II 

© 

(84) 

m  =  (t-l)(J  +  l'  +j  +  l 

(85) 
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For  example,  if  J  =  2  and  1  =  4,  then  there  are  A/  =  4  x  (2+  1)  =  12  match  pointa,  and  the  fifth 
match  point  is  given  by  (x, y)5  =  (i2.pi)  =  (-|  +  2Ax,-|  +  Ay). 


4.S  System  of  Linear  Equation; 

The  systems  of  linear  equations  given  by  Equations  77-  79  may  be  written  in  matrix  form  as 
V*-**  =  Z‘*-’I  (86) 

The  voltage  matrix  Vr'1,i  lsaMxl  vector  whose  elements  are  defined  by 
V*  =  £*,(*,  y)m 

V*  =  £}(x,y)m  m  =  l  ,U  (87) 

VA  = 

The  current  matrix  /  is  a  3  x  1  vector  composed  of  the  unknown  wave  amplitudes  defined  by 
{/n  =  The  impedance  matrix  is  a  M  x 3  matrix  whose  elements  are  defined  by 

=  Crn  <™  [/-  (y,  +  e^»-*>  - 

j(er-l)  \kl  +  ^\JJs  **n  [/”  (y  +  GdziJ 

+  (tr“1)ai//sc"'  cos[/"(!/  +  01  '*  *'  G  iz' ** 

Z*„  =  <VnCOsj^/„  (w  +  I)]  *",*‘*‘ ~ 

{(tf  _1)S^//s  ct"  ,in  [/o  (,/+  l)]  Gd*'W 

+  («r-l)  +  Jjj]  JJs  cos  (y  +  |  Gdz'dl/ 

+ ik&Jc e" -I*  M  e~“"'Gdz'}i',)=l'.„)  (89) 

=  c«neoe^/„^  +  ^je->***‘  + 

|i«eo(cr  -  l)ctn^  y^sin  fl/  +  e-«**'c  dr'dy' 
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-iu«.(£r-l)c,nA//sCos[/»  (i/  +  |)] 

+e.»|-  /  dx') 


t-”'1' G  dz' dj 


(90) 


A  conventional  moment  method  solution  generates  &  square  impedance  matrix,  and  the  unknown 
vector  I  is  found  by  inverting  the  impedance  matrix  Z.  Here,  however,  we  have  a  non-square 
impedance  matrix,  indicating  an  overdetenruned  system  of  equations  (This  is  analogous  to  the 
problem  of  finding  a  ‘'best  fit”  line  given  more  than  two  points  in  a  plane.)  As  suggested  by 
Sarkar  (24),  we  will  solve  for  1  in  a  least-squared-error  sense  by  multiplying  the  matrix  equations 
by  the  conjugate  transpose  of  Z,  Z*T,  yielding 


Z*TV  Z'TZi 


(91) 


Z*TZ  is  a  3  x  3  matrix  which  can  be  inverted  (if  non-singular)  and  I  is  given  by 

i z*Tz\-lz%Tv  =  / 


(92) 


The  “goodness  of  fit”  of  the  solution  is  determined  by  the  residual  r,  given  by 


r=V-ZI 


(93) 


We  may  define  an  error  metric  as  the  average  relative  error,  which  is  given  by 
AvgRdError=ig|^| 


(94) 


4-4  Impedance  Maim  Z**’* 

The  elements  of  Zr,r,J  are  given  by  Equations  88-90.  Unfortunately,  all  of  she  surface  inte- 
gtands  are  singular  at  the  match  point  (x,  y)m  and  must  be  integrated  with  great  care.  Additionally, 
the  integrands  are  functions  of  x  and  y,  which  create  problems  when  we  attempt  to  differentiate 
them. 
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The  solution  is  to  first  define  a  region  Sim  (where  the  subscript  to  indicates  dependence  on 
match  point  location)  about  the  singular  match  point  such  that  the  integrands  in  region  S— 5jm 
are  not  affected  by  the  singularity. 

Then  we  can  exploit  symmetry  and  small-argument  polynomial  expansions  to  analytically 
evaluate  integrals  in  region  Sim  and  numerically  integrate  in  region  5— Sjm. 

Second,  Barksdale  (4)  suggests  the  change  of  variables  u  =  z’—z  and  v  =  j/— y  to  move 
the  singular  match  point  to  the  origin,  allowing  us  to  move  the  z  and  y  dependence  out  of  the 
integrand.  While  the  limits  of  integration  arc  still  functions  of  z  and  y,  this  dependence  is  moved 
from  the  inner  region  to  the  outer  region;  instead  of  integrating  over  regions  S—S\n  and  S\m  in 
xy-coordinaies,  we  now  integrate  over  regions  Sm—S\  and  Si  in  tip-coordinates 

A  more  complete  and  detailed  development  of  impedance  matrices  Zx ,  as  well  as  their 
rather  complicated  form,  is  presented  in  Appendix  B. 


4  5  Two- Dimensional  Scattering  Wtdth 

Once  the  amplitudes  of  the  three  physical  basis  functions  are  found,  we  know  the  equivalent 
polarization  and  conduction  currents  and  can  use  Equations  6$  and  65  to  find  the  bistatic  two- 
dimensional  scattering  width  as  a  function  of  the  observation  angle  phi,  as  depicted  in  Figure  3 

The  form  of  the  current  expressions,  given  by  Equations  72-75,  allows  us  to  evaluate  the 
integrals  of  Equation  65  analytically  (the  evaluations  are  very  simple),  so  that  ow  is  proportional 
to  the  magnitude  squared  of  the  sum  of  three  algebraic  expressions  corresponding  to  the  three 
physical  basis  functions.  This  sum  is 


JTi  ?(*0CO 

( c,„  tin  4  +  ~  **  [cmCi  sin  4  -  c,nC»  cos  4)1 

l  In  -*0,ID  4  J 

where  c  =  (jt0sin4sin(/nh)  -  /„  ax{f„h)}  +  }„ 


28 


<*  =  fj  *•***+  yk0  sin  +  /„  ccs(/„fc)]  -  jk0  sin  <f> 


(96) 


Note  the  tm(x)/z  envelope  which  is  dependent  on  the  strip  width  w  An  advantageous  feature  is 
that  we  can  isolate  the  contribution  of  each  of  the  three  physical  basis  functions  to  the  overall  a^D- 
For  instance,  we  expect  the  scattered  fields  due  to  the  “forced”  wave  to  be  maximal  in  the  specular 
direction. 
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V.  Results  and  Conclusions 


The  equations  developed  in  the  Chapter  IV  were  coded  in  standard  FORTRAN  and  compiled 
and  executed  on  a  micro  VAX  3+  computer  using  the  VMS  version  5  3  operating  system  The 
program,  called  PBFSTRIP,  reads  geometry  and  material  parameters  from  an  input  file,  generates 
the  impedance  and  voltage  matrices,  solves  for  the  unknown  current  vector  (using  Irvin’s  matrix 
inversion  routine  (11)),  and  writes  the  surface  wave  parameters,  solution  errors,  strip  and  dielectric 
currents,  and  bistatic  scattering  width  to  output  files  A  description  and  source-code  listing  of 
PBFSTRIP  may  be  found  in  the  supplement  to  this  thesis. 

In  this  chapter,  we  present  some  results  from  running  PBFSTRIP  and  evaluate  its  perfor¬ 
mance  We  investigate  the  convergence  of  the  solution  as  the  number  and  configuration  of  match 
points  is  varied  We  compare  the  accuracy  of  the  results  and  execution  time  to  data  provided  by  Dr 
Xent  of  Mission  Research  Corporation,  using  a  pulse-Galerkm  moment  method  implementation. 

We  expect  the  physical  basis  functions  to  accurately  represent  the  conduction  and  polariza¬ 
tion  currents  near  the  middle  of  the  scatterer,  but  not  near  the  edges  where  the  scatterer  docs 
not  resemble  an  infinite  structure.  Additionally,  the  physical  basis  functions  do  not  mode1 
diffractions  from  the  corners  of  the  dielectric  slab  and  the  edges  of  the  conducting  strip. 

We  found  that  only  the  Z*  n  impedance  matrix  needed  to  be  included  in  PBFSTRIP  and 
that  Zjbn  and  Z^n  contributed  nothing  to  the  correct  solution.  We  surmise  this  is  because  the 
most  important  boundary  condition  that  must  be  met  is  that  of  vanishing  tangential  electric  field 
on  the  conductor,  £j(x,  -£)  =  0,  which  only  the  Z^n  impedance  matrix  addresses  directiy. 

5.1  The  Run  Matrix 

We  wish  to  observe  the  behavior  of  the  results  as  the  input  parameters  vary.  The  input 
parameters  consist  of  relative  permittivity  crt  slab  width  tv  and  thickness  A,  and  angle  of  incidence 
0,  measured  counterclockwise  from  the  positive  y-axis.  We  thus  generate  a  run  matrix  by  defining 
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Table  1.  Run  Matrix 


a  baseline  configuration  and  perturbing  the  input  parameter*  individually  about  the  baseline  value 
The  run  matrix  is  shown  in  Table  5.1. 


5.2  Convergence  of  the  Solution  for  Various  Match  Point  Confioumltons 

Although  PBFSTRJP  is  written  to  accept  match  points  both  within  the  dielectric  and  on 
the  conductor,  we  find  experimentally  that  only  the  match  points  on  the  conductor  are  significant 
This  is  probably  due  to  the  dominance  of  the  strict  boundary  condition  on  the  strip  (vanishing 
Ej)  The  question  remains  as  to  how  many  match  point?  to  take  on  the  conductor,  and  how  they 
should  be  distributed 

We  ran  several  cases  at  normal  incidence,  0  =  0* ,  a  geometry  that  will  give  us  forward  and 
reverse  surface  waves  (FSW  and  RSW,  respectively)  of  equal  amplitude  To  achieve  this,  we  find 
we  must  have  a  symmetric  match  point  distribution  with  respect  to  the  left  and  right  edges  of  the 
strip.  Ibis  is  equivalent  to  weighting  the  FSW  and  RSW  equally. 

Using  symmetric  distributions,  we  wish  to  determine  the  effects  of  different  numbers  of  match 
points,  and  how  close  to  the  edges  we  may  place  them.  We  repeated  each  case  in  the  run  matrix 
with  5,  9,  17,  and  21  match  points  on  the  conductor  to  test  for  convergence  of  the  solution.  The 
results  for  the  baseline  Case  0  are  presented  in  'Ifcble  5.2;  the  other  cases  show  similar  trends 
"Range”  refers  to  the  portion  of  the  strip  over  which  the  match  points  are  distributed  and  CN  is 
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H 
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Wave  Amplitude 

Avg  Rel  Error 
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H 
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0.001163^100.4® 

19  6 

0.005103 

HI 
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1.0701  —  22  4* 

0  00047U-78  6' 

0.001491Z116.1® 

17.7 

0  012775 

Table  2.  Baseline  Case  Results  for  Various  Match  Point  Configurations 


the  condition  number  of  the  normalized  3x3  impedance  matrix 


Several  conclusions  may  be  drawn  from  Table  5  2.  First,  the  forced  wave  amplitude  is  insensi¬ 
tive  to  match  point  configuration,  but  the  FSW  and  RSW  amplitudes  fluctuate  moderately.  Second, 
the  RSW  amplitude  is  greater  than  the  FSW  amplitude,  which  is  not  surprising  if  we  think  of  the 
incident  wave  “feeding”  the  RSW  at  the  right  edge  of  the  slab  more  strongly  than  the  FSW  at  the 
left  edge.  Third,  the  normalized  impedance  matrix  is  numerically  stable  (low  condition  number) 
for  all  match  point  configuration  Fourth,  the  average  relative  error  (see  Equation  94)  increases 
with  the  number  of  match  points,  but  does  not  become  unacceptably  large  The  overall  conclusion 
is  that  the  physical  basis  function  solution  converges  even  for  relatively  few  match  points 


5.S  Accuracy  of  the  Solution 

Convergence  of  the  solution  is  only  desirable  if  it  converges  to  the  correct  values.  We  use 
data  provided  by  Dr.  Kent  from  a  pulse*Galerkin  code  as  truth  data,  and  look  for  agreement 
both  m  the  currents  and  in  the  scattered  far-fields.  All  cases  defined  in  Table  5  1  were  run  for 
M  =  5,9,17,  and  21,  and  the  currents  and  far-fields  for  the  PBF  and  pulse-Galerkir  solutions 
overlaid  in  a  common  format.  We  compare  PBFSTRIP  and  pulse-Galerkin  results  for  the  extremes 
of  M  =  5  and  M  =  21,  these  results  are  shown  in  Figures  5  through  14.  Each  figure  is  organized 
as  fellows  The  upper  left  graph  shows  the  conduction  current  Je  along  the  conducting  strip.  The 
upper  right  and  lower  left  graphs  show  the  r-directed  and  p-direeted  polarization  currents,  Jff  and 
J/f  respectively,  along  the  middle  of  the  dielectric  slab.  The  lower  right  graph  shows  the  bistatic 
scattering  width  where  the  bistatic  angle  d  is  measured  counterclockwise  from  the  positive 
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|J„’|  (mA/m1)  |J,|  (mA/m) 


Figure  5.  Case  0:  cr  =  5.00,  w  =  5  QAq,  h  =  0  05 A o,  0  =  —30® 

x-axis. 

In  Figure  5  we  see  that  the  conduction  current  Jc  predicted  by  the  PBF  and  pulse-Galerkin 
solutions  grossly  agree,  but  there  are  differences  on  the  order  of  10%  at  five  places  on  the  strip. 
Additionally,  the  pulse-Galerkin  Je  correctly  vanishes  at  the  strip  edges  while  the  PBF  Jc  does 
not.  The  pulse-Galerkin  is  an  almost  constant  10-20%  below  the  PBF  ,  but  the  locations 
of  the  “hills  and  valleys"  agree  well.  Again,  the  pulse-Galerkin  Jf#  disappears  at  the  edges  while 
PBFSTRIP  fails  to  predict  this.  The  pulse-Galerkin  and  PBF  J*f  agree  very  closely.  The 
bistatic  scattering  widths  predicted  by  the  two  solutions  agree  very  closely  in  the  number,  locations, 
and  widths  of  the  main  and  side  lobes.  Both  methods  predict  the  main  lobe  in  the  specular  direction 
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10 


|Jj  vs  X 


15 


|J„’j  vs  X 


Figure  6.  Case  1:  cr  =  2.00,  w  =  5  OAo,  h  =  0  05  Ao,  9  =  —30* 

of  £  =  120*.  The  two  solutions  overlay  almost  exactly  (within  *-2  dBA)  for  <£  between  30°  and 
150*,  and  are  within  3  dB  for  the  near-grazing  bistatic  angles  0*  <  $  <  30°  and  150*  <  ^  <  180® 
Notice  that  for  Case  0,  M  =  5  agrees  more  closely  with  the  pulse-Galerkin  code  than  M  =  21  in 
predicting  the  far-fields. 

Case  1  is  illustrated  in  Figure  6.  Permittivity  cr  is  small,  and  the  polarization  currents  are 
correspondingly  lower  than  in  Case  0.  The  strip  conduction  current  Je  predicted  by  the  PBF 
and  pulse-Galerkin  solutions  show  the  same  rough  equivalence  as  seen  in  Case  0  The  scattering 
widths  predicted  by  the  two  methods  overlay  almost  exactly,  although  the  far-field  PBFSTR1P 
nulls  are  deeper  than  the  pulse-Galerkin  nulls  Otherwise,  the  largest  disagreement  is  about  5  dBA 
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Figure  7  Case  2:  er  =  8.28,  w  ~  5.0Ao,  h  =  0  05Ao,  0  =  —30° 

at  $  =  20° 

Case  2  is  illustrated  in  Figure  7.  Permittivity  cr  is  large,  and  the  polarization  currents  are 
much  higher  than  in  Case  0.  The  strip  and  polarization  currents  predicted  by  the  two  methods 
show  the  same  characteristics  as  noted  for  Case3  0  and  l.  The  far-field  disagreement  is  about  4  dBA 
for  0*  <  0  <  45*  and  150*  <  ^  <  180*  and  the  outermost  nulls  are  displaced  about  4*:  otherwise 
the  fax- folds  overlay  almost  exactly. 

Case  3  is  depicted  in  Figure  8.  Strip  width  w  is  small,  and  the  fur-field  lobes  are  corre¬ 
spondingly  w.der  than  in  Case  0.  The  PBFSTRIP  and  pulse-Galetkin  conduction  and  polamat:on 
currents  show  qualitatively  poorer  agreement  than  for  the  previous  cases,  although  M =21  agrees 
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Figure  8.  Case  3:  c,  =  5.00,  ty  =  3  OAo,  h  =  0.05Ao,  0  =  —30* 

mere  closely  than  M  =  5.  The  far-fields  agree  very  well,  for  M=21  the  maximum  disagreement  is 
less  than  2  dBA  over  the  entire  0°  <  ^  <  180*  bistatic  range.  When  M  =  5,  however,  the  far-fields 
disagree  by  as  much  as  10  dBA  for  0°  <  ^  <  35°. 

Case  4  is  shown  in  Figure  9  Strip  width  w  is  large  and  the  far-field  lobes  are  correspondingly 
narrower  than  in  Case  0.  The  PBFSTRIP  and  pulse-Galerkin  currents  are  roughly  the  same.  The 
far-fields  overlay  almost  exactly  except  near  grazing  (0*  <  6  <  45*,  155°  <  ^  <  180°).  However, 
there  may  be  a  problem  with  the  pulse-Galerkin  solution  in  this  case,  as  evidenced  by  the  very 
deep  cull  at  ^  =:  8*;  all  other  cases  have  a  nearly  constant  scattering  width  when  0°  <  <£  <  30* 
This  null  is  suspect  because,  due  to  the  large  width  of  the  slab,  the  pulse-Galerkin  cell  size  had  to 
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Figure  9  Case  4:  cr  =  5.00,  ty  =  IO.OAq,  h  ~  0.05Ao,  0  ~  —30° 
be  raised  above  the  convergence  threshold  of  A*/10  to  A*/9. 

Case  5  is  illustrated  in  Figure  10.  The  slab  is  half  as  thick  as  in  Case  0.  As  before,  the 
PBFSTRIP  and  pulse-Galerkin  currents  are  roughly  equivalent,  although  the  PBFSTRIP  J*#  is 
nearly  twice  as  large  as  the  pulse-Galerkin  Jff.  The  far-fields  overlay  almost  exactly  except  for 
0°  <  ^  <  50®,  where  the  PBFSTRIP  scattering  width  is  as  much  as  5  dBA  lower  than  the  pulse- 
Galerkin  scattering  width.  Also,  the  PBFSTRIP  nulls  are  10-15  dBA  below  the  pulse-Galerkin 
nulls. 

Case  6  is  shown  in  Figure  11.  The  dielectric  slab  is  twice  as  thick  as  in  Case  0,  although  not 
thick  enough  to  allow  higher-order  surface  wave  modes  to  cut  on  For  the  pulse- 
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Figure  10.  Case  5:  er  —  5.00,  u>  =  5.0Ao,  h  =:  0.025Aq,  $  =  ~30® 

Galerkin  solution,  we  are  forced  into  using  a  cell  size  of  A^/5,  so  the  pulse-Galerkin  results  are 
suspect  Still,  the  currents  predicted  by  the  two  methods  agree  fairly  well,  although  the  currents 
appear  to  be  “shifted”  spatially  on  the  x-axis  by  about  0.25Aq  The  far-fields  agree  fairly  closely, 
except  in  the  region  25*  <  4>  <  50°  where  the  PBFSTRIP  scattering  width  is  10  to  15  dB  below 
the  pulse-Galerkin  prediction 

Case  7  is  depicted  m  Figure  12.  The  slab  is  illuminated  at  normal  incidence  ($  =  0*),  and 
the  main  lobe  is  in  the  backscatter  =  90*)  direction.  Both  PBFSTRIP  and  the  pulse-Galerkin 
solution  show  the  expected  left- to- right  symmetry  in  the  currents.  This  geometry  allows  us  to 
clearly  see  the  effects  of  the  three  physical  basis  functions;  the  forced  wave  currents  are  exclusively 
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Figure  11.  Case  6:  cr  =  5.00,  w  =  5.0Ao,  h  =  0.10Ao,  0  =  -30° 

x-directed  and  their  phase  is  constant  with  respect  to  x,  and  the  surface  waves  form  a  pure  sinusoid 
pattern  “riding”  on  the  forced  wave. 

Another  current  component  is  evident  in  the  pulse-Galerkin  Jc,  probably  due  to  diffraction 
effects  from  the  strip  edges.  The  far-fields  agree  for  45*  <  ^  <  135°  but  differ  significantly  for 
bistatic  angles  closer  to  grazing.  PBFSTRIP  predicts  two  nulls  in  the  angular  region  0*  <  4>  <  45° 
while  the  pulse-Galerkin  solution  predicts  -5  ±  2  dBA  for  0 *  <$<  40®. 

Case  8  is  shown  in  Figure  13.  The  slab  is  illuminated  at  60*  from  normal  ( 6  =  -60°)  The 
PBFSTRIP  and  pulse-Galerkin  methods  both  predict  a  current  distribution  consisting  of  a  low 
spatial  frequency  component  with  a  “ripple”  riding  on  it.  This  can  be  explained  by  a  strong  reverse 
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Figure  12.  Case  7:  cr  =  5.00,  w  =  5.0Ac,  h  =  0.05Ao,  0  =  0° 

surface  wave  “beating  against”  the  forced  wave,  and  the  iow  spatial  frequency  is  determined  by  the 
difference  in  their  r-directed  phase  constants.  This  corresponds  to  a  “beat”  wavelength  of  5  76Ao, 
corresponding  fairly  closely  to  the  peak- to- valley  distance  of  about  1.7Ao  shown  m  Figure  13  The 
far- fields  agree  fairly  well  except  for  0*  <  ^  <  45*  where  the  two  solutions  disagree  by  as  much  as 
10  dBA  for  0*  <  *  <  35*. 

Case  9  is  illustrated  in  Figure  14,  The  slab  is  illuminated  at  very  near  grazing  incidence, 
0  =  —85*.  We  see  the  same  “beat"  phenomenon  as  in  Case  8,  but  now  the  “beat"  wavelength 
is  23.3Ao,  much  longer  than  the  strip  width  The  PBFSTRIP  and  pulse- Galerkin  currents  are 
roughly  equivalent  except  for  a  large,  but  uniform,  disagreement  in  the  The  two  scattering 
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Figure  13.  Case  8.  cr  ~  5.00,  w  =  5.0Aq,  h  =  0.05Aq,  0  =  —60° 


width  predictions  show  the  same  trends  and  lobing  patterns,  but  the  PBFSTRIP  scattering  width 
prediction  is  1  to  10  dB  below  the  pulse-Galerkin  prediction. 


S  4  Execution  Times 

The  real  advantage  of  using  a  physical  basis  function  solution  is  the  reduced  computer  time 
and  memory  required  when  compared  to  subsectional  solutions,  such  as  the  pulse-Galerkin  tech¬ 
nique  considered  here.  The  execution  times  for  the  ten  test  cases  given  in  Table  5.1,  as  well  as  the 
dimension  N  of  the  pulse-Galerkin  impedance  matrix,  are  shown  in  Table  5  4  Both  codes  were 
run  on  a  microVAX  3+. 
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Figure  14  Case  9:  er  =  5.00,  w  =  5.0Ao,  h  =  0  05A0,  $  =  -85* 

The  PBF  solution  executes  about  10  times  faster  than  the  pulse»Galerkm  solutions  More 
importantly,  the  PBF  solution  does  not  depend  strongly  on  the  size  of  the  scatterer,  while  the 
pulse-G&lerkin  solution  execution  time  increases  greatly  with  scatterer  size. 

The  pulse-Galerkin  solution  uses  most  its  execution  time  in  solving  the  NxN  impedance 
matrix  by  Gaussian  elimination  As  reported  by  Woodworth  (27),  this  time  is  proportional  to 
AT3.  The  pulse  size  should  be  no  larger  than  A^/10,  where  Xj  =  A0/v/c7  is  the  wavelength  in  the 
dielectric,  so  N  is  proportional  to  the  size  of  the  scatterer.  Additionally,  “heavy”  dielectrics  (large 
c-)  require  more  pulses  than  light  dielectrics,  so  N  also  increases  with  increasing  cr. 

The  PBF  solution  expends  most  of  its  execution  time  in  performing  numerical  integrations 
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Case 

PBF 

CPU  (sec) 

Pulse*Galerkin  1 

u 

K22XCS] 

mm 

23.5 

kMii 

ESI 

■  mm 

, 

83  8 

25.0 

1617.5 

5 

23.7 

300 

285  1 

6 

24.3 

420 

799.8 

21.2 

300 

293.8 

23.9 

300 

294.6 

23.8 

300 

296.4 

Table  3.  Run  Times:  Physical  Basis  Function  vs  Pulse-Galerkin  Solutions 


to  fill  the  3  x  M  impedance  matrix.  Negligib'e  time  is  spent  ‘‘normalizing”  the  impedance  matrix 
(multiplying  by  its  conjugate  transpose)  and  solving  the  resultant  3x3  system  of  equations.  The 
integration  time  increases  with  scatterer  width,  which  is  evident  in  Case  3  but  not  in  Case  4,  where 
the  PBFSTRIP’s  integration  algorithm  increased  the  Gaussian  quadrature  node  spacing  rather 
than  the  number  of  nodes. 


5  5  Components  of  ike  Far-Fttld 

The  use  of  physical  basis  functions  allows  us  to  gain  some  insight  about  how  the  different 
scattering  mechanisms  contribute  to  the  total  scattering  width  as  a  function  of  bistatic  angle  <ff 
Figure  15  shows  the  total  scattering  width  a*  ‘  the  contributions  from  the  forced  wave,  FSW,  and 
RS W  for  the  baseline  Case  0,  M  =  5  We  see  that  the  forced  wave  strongly  dominates  the  scattering 
width  for  35°  <  <t>  <  150*,  being  nearly  30  dB  above  the  surface  wave  contributions  at  the  main 
lobe.  The  surface  wave  contributions  dominate  the  scattering  width  at  near-grazmg  angles;  the 
FSW  dominates  for  0®  <  4>  <  35°  while  the  RSW  dominates  for  150®  <  4>  <  18G®.  Another  effect  of 
the  surface  waves  is  to  counter  the  deep  forced  wave  nulls  at  25®,  46®,  59®,  and  154®. 
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PBF  Components  of  the  Bistctic  Scottering  Width 


Figure  15.  Baseline  Case,  5  Match  Points:  cr  =  5.00,  xv  =  5.0Ao,  h  =  0.05Ao,  0  =  —30° 

5.6  Conclusions 

The  physical  basis  functions  provide  a  reasonably  complete  and  very  efficient  set  in  which  to 
expand  the  currents  and  predict  the  far-fields  for  the  dielectric-coated  conducting  strip 

The  PBF  solution  is  relatively  insensitive  to  the  number  of  match  points  M  so  long  as  M  >3 
Match  points  on  the  conductor  work  best,  match  points  inside  the  dielectric  slab  do  not  contribute 
significantly  to  the  solution  A  set  of  match  points  chosen  symmetrically  gives  the  best  results 

The  PBF  solution  yields  only  roughly  correct  currents  on  the  conductor  and  in  the  dielectric 
slab,  but  gives  very  accurate  far-field  predictions  for  incident  and  bistatic  angles  away  from  grazing 
At  grazing  bistatic  angles,  diffractions  from  the  slab  corners  and  strip  edges  become  significant, 
and  these  contributions  are  ne  t  contained  in  the  physical  basis  functions 

The  PBF  solution  runs  an  order  of  magnitude  faster  than  a  standard  pulse-Galerkin  moment 
method  solution  for  scatterers  cf  even  modest  dimensions  (5A<>  wide  x  0  05Ao  thick).  The  PBF 
solution  can  be  applied  to  large  scatterers  where  the  pulse-Galerkin  solution  simply  cannot  be  used 
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If  used  with  care,  the  speed  of  the  PBF  solution  could  allow  us  to  find  the  far-fields  for  a  range  of 
different  frequencies  and  apply  Fourier  analysis  to  predict  a  time-domain  response  of  the  scatterer 

5.7  RecommendaUons 

While  the  PBF  solution  and  software  developed  in  this  thesis  effort  predicted  the  two- 
dimensional  scattering  width  of  a  dielectric-coated  conducting  strip  fairly  well,  it  did  not  work 
well  when  the  incident  or  bistatic  angle  was  near  grazing  Unfortunately,  agreement  with  the 
pulse-Galerkin  currents  was  less  than  ideal  In  short,  it  was  not  perfect  and  can  certainly  be  im¬ 
proved  upon;  additionally,  its  scope  could  be  extended.  Several  aieas  where  future  research  could 
explore  follow 

First,  PBFSTRIP  could  be  extended  to  include  lossy  and/or  magnetic  material.  Lossy  ma¬ 
terial  would  be  relatively  easy  to  incorporate;  the  only  difficulty  would  be  finding  the  surface  wave 
phase  constants  from  the  transcendental  equation  (Equation  124)  when  cr  is  complex  Muller’s 
algorithm  (16),  as  implemented  by  the  IMSL  subroutine  ZANLY,  could  be  used  for  this  task  The 
inclusion  of  a  magnetic  material  (p<*  po)  would  be  slightly  more  difficult  to  implement  as  the 
integral  equations  would  then  include  an  additional  surface  integration  of  the  equivalent  magnetic 
polarization  current  ,  as  given  by  Equations  46,  51,  and  55. 

Second,  ‘.he  boundary  conditions  at  the  vertical  edges  of  the  dielectric  slab  and  conducting 
strip  are  not  satisfied;  remember  that  the  physical  basis  functions  come  from  infinite  structures 
that  have  no  vertical  boundaries  The  results  clearly  show  that  PBFSTRIP  does  not  predict  the 
necessary  vanishing  x-  directed  currents  at  the  vertical  boundaries.  It  would  certainly  be  valid 
to  devise  a  hybrid  solution  in  which  physical  bases  are  used  in  the  middle  of  the  structure  and 
a  few  rationally  chosen  subsections!  basis  functions  are  used  at  the  edges.  An  example  could  be 
one-dimensional  piecewise  sinusoids  on  the  strip  and/or  two-dimensional  rooftop  functions  near  the 
edges  of  the  dielectric  slab.  Dr.  Kent  used  this  technique  effectively  for  the  similar  problem  of  a 
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dielectric  slab  on  an  infinite  ground  plane  illuminated  by  a  TM*  plane  wave. 

Third,  physical  models  of  a  dielectric-coated  conducting  strip  could  be  fabricated  and  its 
bistatie  scattering  width  could  be  measured  in  an  anechoic  chamber.  Since  a  true  two-dimensional 
model  obviously  is  not  physically  realizable,  a  three-dimensional  strip  of  length  L  could  be  made 
and  its  three-  dimensional  radar  cross  section  <730  measured  Then  the  two-dimensional  scattering 
width  <720  is  given  by  (3.578) 

aw  ~  £  >  *0  (97) 
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Appendix  A.  Physical  Basis  Functions 


The  physical  basis  functions  are  those  wave  components  which  exist  on  a  dielectric-coated 
conducting  strip  of  infinite  width.  In  the  limit  as  tu-*oo,  the  two-dimensional  scattering  geometry  of 
Figure  1  degenerates  into  a  one-dimensional  dielectric-ewered  ground  plane,  as  shown  in  Figure  16. 


A.J  Forced  Wave 

Consider  the  geometry  formed  by  a  ground  plane  covered  with  a  dielectric  slab  of  height  h 
Let  a  TE 2  uniform  plane  wave  be  incident  or*  the  dielectric,  as  shown  in  Figure  17,  which  sets  up 
two  uniform  plane  waves  within  the  dielectric  and  a  reflected  uniform  plane  wave 

The  elec*nc  fields  in  free  space  are  composed  of  incident  (E/)  and  reflected  (Er)  components 
The  fields  in  the  dielectric  arc  composed  of  “downward”  (E*1)  and  “upward”  (E*3)  components. 
By  Snell’s  Laws  (3191),  0  =  0r*  and  kg  sin  $d  =  k0nn6  or  y/c? sui  64  r-  sinfl.  We  choose 

a  reference  point  at  (0,  £)  and  arbitrarily  force  the  fields  to  have  zero  phase  at  that  point.  We  then 
write  the  electric  fields  explicitly  as 

E*(z,y)  =  Coe*4*tin0)  (98) 

E '(*.»)  =  Xe-iM*"”,+<»-*>c“,)(S*w0-fain0)  (99) 
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Figure  17.  “Forced”  Wave  in  Dielectric-Covered  Ground  Plane 


EJ1(x  ,y)  =  Be--'*'<*'m''-<'‘I>e“,')(xcos^  +  ysin9j)  (100) 

E  "(r,y)  =  Ce-'*‘(*““*<+(»-I>c“'')(xcosflJ-y«in^)  (101) 


Without  loss  of  generality,  we  let  E0  =  1. 

The  magnetic  fields  are  related  to  the  electric  fields  through  the  impedance  of  the  appropriate 
medium,  and  are  given  by 


H‘(*.») 

*7# 

(102) 

Hr(x,y) 

Vo 

(103) 

H'^x.y) 

Vi 

(104) 

H"(x,y) 

_  #<+(»- 
Vi 

(105) 

where  Vo  =  \/p*/c0  and  w  =  \jiulu  =  We  we  interested  in  the  fields  withm  the 

dielectric,  which  are  given  by 

(106) 

H'(*,y)  =  H*(*»y)  +  H«(z,y) 

We  now  apply  the  boundary  conditions  and  solve  for  the  unknown  constants  A,  B,  and  C  The 
first  boundary  condition  requires  that  the  total  tangential  electric  field  vanish  along  the  dielec- 
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tnc/conductor  boundary.  We  arbitrarily  choose  the  point  (0,—^)  and  aee  that 


Ei'  (o.-^  +  Bf2  coa$d  +  Ct*>t‘ha,,‘ 

or  equivalently 


B4-Ce+,2k<'*co-#<  =0 


(107) 


The  second  boundary  condition  requires  that  the  total  tangential  electric  and  magnetic  fields  be 
continuous  across  the  free  space/dielectric  boundary.  We  arbitrarily  choose  the  point  (0, 4)  and 
see  that 


£*(°’i)+£*(o,l)  =  ^(°.|)+^(°.|) 

m  (°.|) +«:  (°.|)  =  «f( 0.5)+*."  (o.f) 

Plugging  into  Equations  98-105,  we  can  rewrite  this  as 


cos0  +  i4cos0 


B  coedd  +  Ccc»$d 


±__±  -  £_£ 
Vo  Vo  Vi  Vi 


or 


cosQ  cos  6 


A  +  B— -C— 

Vi  Vi 


-1 

1 


(108) 

(103) 


We  put  Equations  107,  108,  and  109  into  matrix  form 


0 

1 

yl 

0 

1 

_$2 l2i. 

CO*  i 

-Sfr 

B 

= 

-1 

1 

a*. 

_2*. 

a* 

C 

1 

(110) 
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Ftom  Equation  106,  we  only  need  to  find  B  and  C.  Solving  Equation  110,  letting  £  =  kfh  cosfy, 
and  noting  that  r)0/i)d  —  y/cri  we  find  the  following  expressions  for  the  B  and  C 


cosdcM 

y/e^co&Ocost  +  jcos0rfsm£ 
— cosd 

^7 cos 0  cos£  +  ;  cos$d  sin ( 


an) 

(112) 


We  can  now  find  expressions  for  the  total  electric  field  within  the  dielectric  using  Equations  100, 
101,  and  106.  First,  let  us  consider  the  x-component.  Taking  advantage  of  cancellations  and 
remembering  Snell’s  Law  of  Refraction,  we  find 


Ei{x,y)  =  + 

—  L}My+$)c<»#4  _c-fM*+$)co«i4l _ eo<g  cm  $4 _ >n# 

l  J  v/e7cos0cos(  +  ,7coe0(f8in£ 


=  j  cos  $d 


2  cos$ 


yftZ  cos  6  cos  £  +  j  cos $i  $in( 


sin  I  kd  cos  $d 


(•*$)]■ 


tx 


(113) 


Similarly,  we  find  the  y-component  of  the  electric  field  to  be  given  by 

=  .  — »«■« -co.  [W,  (y+  «-»*.-»«  (1M) 

^/F7cos0cos£  +  .;co60*8in£  [  \  2/J 


The  z-directed  toted  magnetic  field  in  the  dielectric  may  be  found  using  Equations  104, 105,  and  106 
and  is  given  by 


Hf(z.y) 


iLc-zM*«n #4 -(»-!)»•#*)  __  £e-iM*«n  *<+(»- $)«»♦<) 
1)1  1)  i 


y/cj  _ 2  COS# _ 

*7#  y/fi?cos0cosZ  +  jcos$d*inZ 


(115) 
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A.2  Surface  Waves 


Consider  the  waveguide  formed  by  a  ground  plane  covered  with  a  dielectric  slab  of  height 
h,  as  shown  in  Figure  16.  By  applying  the  appropriate  boundary  conditions,  we  find  that  TU‘ 
(odd)  and  TE‘  (even)  modes  will  propagate,  but  that  the  TAfiJ  mode  is  dominant  and  has  a  cutoff 
frequency  of  sero  For  thin  dielectric  slabs  only  the  dominant  TM{  mode  will 

propagate  unattenuated,  and  the  electric  fields  within  the  dielectric  are  given  by  (8:163-184) 


£«(*»*)  = 

4£sin[/(w+i)]e-^ 

(116) 

E  (x  - 

-4Iico.fr  fu+  4\L-r,s 

(117) 

««*  cosf' 

£«(*,»)  = 

0 

(118) 

By  Faraday’s  Law,  H  =  x  E,  the  TM£  magnetic  field  is  given  by 


H,  =  H,  s  0  (119) 

H,  =  ~Af  cos  |/  (120) 

The  propagation  constants  /  and  g  may  be  found  using  (3:164) 

(/  h)  tan(/  h)  =  ahcr  (l2i) 

f7  +  =  *2  =  fr  k]  (122) 

-a5  +  g7  =  h]  =  u7it,c„  (123) 

Equations  121-123  may  be  combined  into 

(/  h)  tan(/ h)  =  Crhs/(Cr  -  l)tj  -  p  (124) 


Equation  124  is  a  transcendental  equation  in  the  phase  constant  /  whose  solution(s)  must  be  found 
numerically,  although  Richmond  gives  a  graphical  technique  for  real  c,  that  is  accurate  to  about 


three  figures  (18).  /  is  defined  to  be  positive,  but  if  /  is  negative  the  fields  are  unchanged  except 
for  a  sign  reversal  of  wave  amplitude  A.  y  is  defined  to  be  positive  for  a  surface  wave  propagating 
sn  the  +r-direction. 

Furthermore,  if  h  <  then  only  one  solution  for  /  exists  and  lies  within  the  range 

0  <  /  <  -  Once  the  phase  constant  /  is  found,  phase  constant  g  may  be  found  using  Equation  122 

and  the  fields  may  be  found  using  Equations  116-120. 

The  reverse  surface  wave  has  the  same  form  as  the  forward  surface  wave.  Phase  constant  g 
is  negative,  however,  while  phase  constant  /  remains  positive. 

A.S  Phase  Constants  f  and  g 

The  surface  waves  and  forced  wave  have  the  same  form  except  for  the  amplitudes  of  the 
individual  field  components  and  phase  constants  in  the  z  and  y  directions.  We  may  compactly 
write  expressions  for  the  three  types  of  waves  as 


Err.(z,y) 

—  An  Cm  j/i 

’M) 

I  «-'»•* 

(125) 

E,n(z,y) 

—  An  Cjrn  COS  |/ 

"(y+!) 

(126) 

—  An  Cjn  COS  |/ 

"(y+8 

(127) 

where  n  =  1  represents  the  forced  wave,  n  =  2  represents  the  forward  surface  wave,  and  n  =  3 
represents  the  reverse  surface  wave.  The  unknown  constant  An  is  the  amplitude  of  the  nth  wave 
Phase  constants  /„  and  gn  are  given  by  Equation  113  for  n  =  1,  and  by  Equations  122  and  124  for 
n  =  2  and  3. 

The  known  constants  ctn,  c,n,  and  e,„  define  the  relative  amplitudes  of  the  Ett  and  Ht 
components  of  the  n1*  wave  and  are  defined  by  Equations  116, 117, 120, 113, 114,  and  115.  Written 
in  convenient  matrix  form,  they  are  given  by 


C.l 

Cx7 

C,3 

j  COS  04 

Jb- 

)wt4 

i**4 

C,2 

C,3 

— 

sin  04 

Ill 

Ut4 

zlx 

Wt4 

Cii 

Cj  2 

yferh h 

-1 

-1 

(128) 


where  04  —  arcsinfcinfl/^/cIT).  Maxwell’s  equations  (specifically,  Ampere’s  and  Gauss’  Laws)  allow 
us  to  relate  the  field  amplitudes  by 


cxnfn  +  jCynQn  —  =  jkoTfyCt> 

Cyn/n  +  JCznQn  =  0 


(129) 

(130) 
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Appendix  B.  Calculation  of  Impedance  Matrix  Elements 


To  fill  the  impedance  matrix  we  must  evaluate  expressions  involving  the  derivatives  of  surface 
and  line  integrals.  We  will  find  that  the  results  are  combinations  of  surface  integrals  over  the 
dielectric  slab  cross-section,  line  integrals  along  the  edges  of  the  slab  and  along  the  conductor,  and 
integrand  evaluations  at  the  comers  of  the  slab  and  edges  of  the  strip. 


B.J  Impedance  Matrix  Element  Z^n 
Zmn  *3  6*vcn  by  (see  Equation  88) 

»in({+) 


We  introduce  the  change  of  variables  tx  =  x1  -  z  and  v  =  ]f  -  y  and  re-write  Equation  131  as 
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cosK+)eos  {/„») 
-sin({+)sin(/„») 


e-ll.('.+-)J{m  +  tlj)  </ud„ 

^  1*“ +  S]  L  (*V*+<*  +  5>*)  d“  035) 


The  regions  of  integration  are  depicted  m  Figure  18.  The  tilde  on  the  regions  of  integration  indicate 
integration  over  uv-space  and  the  m  subscript  indicates  dependence  on  the  match  point,  (z, ,  y} ) 
Region  S\  is  a  disk  of  radius  R  and  is  constant  with  respect  to  z,  and  y; . 


Zr  _ 

mn 


c,n  sin(^+)e■,'*I• 

-c-jt  [c*n  (*♦ + h) 8in({+)  +  c”i^coeK+)] 


*»(«».•) 


Us  S  CO8^/“,,)€’W*U^02,  (*•>/«*  +  V2)  durfv 

-^r  [c*"  (*• + £)  cos«+>  - 


X 


/»<•.») 


ff  sin + »J) 

■»  JSm-Sl  '  ’ 

_£4T  I'”  (**  +  S)  ,ia!£+)  +  <*aS»C"(W] 
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'  '  '*~  '  •N 

h + £]  4'”“‘  L  (Wus+te+^) du  (i36) 

/finj  denotes  the  surface  integral  over  the  circular,  or  semicircular  if  y;  =  —  region  Si  that 
surrounds  the  singular  origin  in  uv-space.  Iftn/  is  derived  in  detail  in  Appendix  C. 

When  we  evaluate  the  derivatives  of  the  functions  outside  the  integrals  we  find  we  can  write 
Zmn  =  sinfo) 

sinK+)  [C™  «  -  ff»2)  + 


-  gj)  +  i<Vn/nSn]  [sm ((+)//  F|(u,  v)  du  dv  +  CO «({+)//  Fj(u, v)  du  *»] 
-[<V»/»  +  j2cIn}„]  [»in(?+)&//Fi(u,v)  dudv  +  cos(£+)£  Jf  f2(u,  v)  dudv] 

'  [cos«+)^//F,(«1v)</ti</ti-8in({+)i//F3(ti,v)dudK] 

+c,„  [«*(£+)jfjj/./'F,(u,v)  dudv  -  sinK+)jfjj  Jf  F7(u,  v)  dudv ] 

+c„  [sin({+)j^r //  Fi(u,v)  dudv  +  coe((+)£,  ff  F2(u,v)  dudv] 

/jsWd„  (.37) 

In  Equation  137,  functions  F\  and  F?  are  integrated  over  region  Sm—Si.  We  evaluate  the  derivatives 
of  the  integrals  using  a  simplified  form  of  Leibnitz’  Rule,  applicable  when  the  integrand  is  not  a 
function  of  x  (1:478) 

ij^FiOdo  =  F[j(i)]^S(r)-F[/(*))^/(r) 
and  find  we  can  write 

ill  F,(u,v)dudv  =  (138) 

&JJ  F}(u,v)dudv  =  -  c_i,‘ ?/„„/]  (139) 
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h!z^)iu  = 

Ty!Lj^'v)iudv  = 

HL-l'{u'v)dudv  = 

—  [  f  Fi(ti,v)  dudv  = 

&IL^d^  = 

= 

TT'!Lj:(u'v)dudv  = 

JjL  j6F3{«,y)du  = 

HLj:iu'v)dudv  = 


*"-*■  [«J«.?//<2)(t0(++))  -  e-».?/,W(to(-+))]  (140) 


«*(£+) 


-cos((-) 


sin({+) 


+«in((_) 


C06({+)/m„i  -  CC*K_)/mn| 

-  «m(4+)7m„k  -  sin(f-)/m„i 

<J»*?w':)(to(++)) 

>(*<,(++» 

-e-»-*//£2,(*o(-+)) 
<J»-*ff<2)(to(+-)) 


-e”‘*-k0 


G+») 


(++) 


t*'' 


'"'9  {iSnlmnc  -  io(f  + 

-‘-“•9  -  *.)7™k] 

c‘ ** 9  [jgnImn3  —  io(|  +  ^iKmm) 
-c~»"9  [jg„lmnJ  +  i0(|  -  ii)/m«,| 

;s„42>(M++)) 

■  /«,  .  >»i”(W++ 

-'-ft?)— 

;jn4,)(*«(-+)) 

[s>n(£+)/„,„k  +  sin(£-)/m„r] 

—to  {c06(^)/mnm  4"  COd(^_)/nmo) 


e»' 


c>>‘9 


-t-il-9 


(141) 

(142) 

(143) 

(141) 


}  (1«) 


(1«) 


(147) 


(148) 


(149) 


(150) 
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4-Z  If  F2\u,v)dudv  =  -fn  [coe({+)Jmn*  -coe(£_ 

or  J  Js^-Sx 

+^0  [sin(^+)/(nnm  “  WD({-)/mnn] 


(151) 


where 


Anna  =  /  /  Fj(tl,  v)  du  dv 
Ann*  =  //  Fi(utv)dudv 

Anne  =  J  cos (/„u)tf<2)  (i0^( |  +  r,)2  +  *>2)  <»« 

/m»i  =  J  sin *> 

Imn.  =  ”  COS (/„»)«?>  (io^f-^  +  O2)  do 

Zmr./  =  J  ^  sm(/„j)//*J)  -l,)2  +  o2)  do 

/mn* =  U-„cos{Uv)  vtr+'sF+s* 

=  j  ^  cos{fnVy  -y^—.  -do 

ft-”  ■ ,,  , j/P’ (*ov/(?  +  *.)4  +  »4>  , 

=  i.w,8(/)  * 

/*-"  .  „  ,«{”  (*«>/(? -*i)s  +  »’)  , 

7m„,  =  do 

Un  =  /*’*'  C-»-“^2)  du 

Imn,  =  f’1'  (to^’  +  fl  -»)*j  du 

*^*0 


J  +  v, 


'5  +  (j +  »)’  j  * 


) 


/?-*«  , 

Annm  —  /  ®  ^  / - 

y  “2  +  (2  +  %)2 

=  /*“'  - ffj2)  (koJu’  +  £-Vi)A 

•'-?-*•  V“s  +  ^-V))’  V  V  / 

(++)  =  '/(f  +  *.)’  +  (3  +  w)2 


<fu 


(152) 

(153) 

(154) 

(155) 

(156) 

(157) 

(158) 

(159) 

(160) 
(161) 
(162) 

(163) 

(164) 

(165) 

(166) 
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(167) 


(+-)  =  \f(  f +  *.)’  + (5 -ft)5 

(-+)  =  \J(  +  ^  + 

(— )  =  )/( f -*.)*+ (3 -»)> 


(168) 

(169) 


Surface  integrals  7mna  and  /mn&  and  line  integrals  /mne  through  Imnn  may  be  evaluated  numerically 
except  that  when  the  match  point  lies  on  the  conductor,  line  integral  Imnk  must  be  integrated 
analytically  in  the  interval  JuJ  <  R  and  /mnm  must  be  evaluated  in  the  limit  as  y; 

Employing  the  same  methods  as  used  to  calculate  I*’**  in  Appendix  C  and  requiring  R  < 
we  may  write 


/—R  >/t 

Jr  ‘'”‘“wol)(*o«)<f«4  2  J  co6 {gnu)H¥\keu)du 

/—  R  /?“*• 

4/ 

2  /  4/  *V"**«JJ)(*o«.)*i 

J-f-z.  Jr 


S  S2i+2t+i 

{a*  -  jPt  fln(t0H)  “  /racI2»  +  2fc  +  1]} 


(170) 


Evaluating  Imnm  when  the  match  point  lies  on  the  conductor  must  be  done  in  the  limit  as 
|  We  find  the  only  contribution  comes  from  the  neighborhood  of  the  singular  point  and  is 
given  by 

2 

Um  Imnm  “jV*  (171) 

i, —  *  *o 

When  we  plug  Equations  138-151  into  Equation  137,  we  find  that  Z*An  is  given  by 


35m  =  Czn  sin(*f)e-"-,‘ 
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-^r-e~iur'c,„kl  [«»(4+)imn.  +  sin ((+)/mnt  +  //,„,] 

+  “  C  ^*r'Cyni?n  (/mni  COs(/n/l)/fliiij] 

gr — 1  |  p'-* 45)(M++))  - Hf >(to(-+))] 

4;  }  - coe(fnh)  [^f-f  4J)(io(+-))  -  e-».»4’>(io(— ))] 

cr-l  ,  f  «n«+) [*<»•»(? +  x.)/™,  +  c-^-»(? -*.)/».»] 

+  -T^-Cx„fc0  < 

3  [  +casK+)[^t(?  +  x1)/fnm+e^-?(  *-*.)/««,] 

+£~2  [^'-?//<5)(m++))  -  «--,,-?4J>(Jto(-+))] 


cin*70 

4 


■rffcf  +  ,  )"P(*o(++»  .  _  x  )^(toi-+» 

*  (2+I,)  (++)  +e  (2  x,)  (-- r) 


(172) 


5.2  Impudence  Matrix  Element  Z^n 

Impedance  matrix  element  Z*,n  is  derived  from  Equation  78  and,  after  implementing  the 
same  change  of  variables  as  in  the  previous  section,  may  be  written  as 


Zln  ~ 


CjnCOsK+Je'"-*1 
(Cr  ~  l)Crn  d1 


JL±+JL 

e~M.(«.+«)/?(»)  (k9\/ tiS  +  t>»J  dudv 

1  (,.^575)  *,* 

»in({+)»in(/„v) 

J6  c->»-<*-+*>4,)  ^>‘,  +  to  +  |),j  du 


4  j  dxdy . 

»in({+)cas(/„v) 
+  cos({+)»in(/„c) 


4^^  5x#y 


(173) 


CO 


Grouping  terms  and  evaluating  the  partial  derivatives  outside  the  integrals  yields 


=  <v„  coe({+) 

4j~e  !l'r  (“V"  (*»  —  f"2)  ~Jc*nfnSn] 

i'»n(to-/«,)-ic*n/nJn)  [«»(£+)// Fi(u,  u)  dudtl  -  SinK+)  //  Fj(tl,v)dudu] 
+Crn/n  [«*(6t)£  //  (“>■')  dudo - *">({+)£  //  F2(u,  o)  dudo] 

-  pc,n/„  +  jc,„g„]  [sin({+)d|;  //  F,(u,  v)  du do  +  cosfo.)^  //  F2(u,  o)  du do] 

+c„  o)  du  dv  +  co*Kx)  //  F2(u,  o)  du  du] 

4c,„  [cos({+)^//Fi(u,»)dudu-sin({+)^/j  Fj(u,u)duduj 

+^e'j,"‘  [-*•£+ &]  L  Fa[u'v,)  du  (m) 

Plugging  in  the  expressions  for  the  partial  derivatives  of  the  integrals,  as  developed  in  the  previous 
section,  we  find 


Z*mn  =  c*n  cos(e+) €-"•*• 

“  c  ""'Cinkl  [cosfo.)  7mn»  8in(£fr  )  Imna  +  ^tng] 

Cr-l_  I  «os((+)  | 

“  ^  ®  *  fcjrn/n  "f  Jtfxnffn]  8jn(/fjh)/m«I 

+  ^c„sin(/„l>)  [e>f*?43)(i0(+-))  _  e-»-f  ff<5>(i0(— ))] 

.  (gr  ~  l)Cyn^O  ~  cini?n>? 0  * 

t  . .  c  itnnm 

COS  (fntymr.i 

C,t,<b  fh  \  [  i,.4^j8)(t»(++))  -(t.T  — h)) 

4  V*  *V[  (++)  (-+) 


(175) 
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B.S  Impedance  Mainz  Element  Z^n 

The  impedance  matrix  elements  associated  with  the  z  integral  equation  may  be  derived  from 
Equation  79  and  are  given  by 


z-m„  =  c,„cos({+)e-"-*' 

uco(cr  -  pc,„  s  rr  rr 
4  8y  JJs.-S,  JJs, 

sin(£+)cos(/r,u) 

+  cos(f+)sin(/„») 


(WuJ  +  u*J  dudv 


UI«o(Cr  -  l)c,„  S 

e-».(r.+-)W(J)  (toV/UJ  +  U*}  dudu 


4  ax 

cos(4+)co «!/,«) 
-smK+)sin(/»c) 


(176) 


Again,  evaluating  the  derivatives  outside  the  integrals,  we  find 


c.n  cos({+;e_/»-*' 

4i?o 

+  jc,r>9n]  [«»({,  ) //  Fi(u,  u)  du  dv  -  sm((+)  ff  F2(u, u)  du du] 
-<V»  [«»(?+)&  //  Fi(u,  v)  du  du  -  *in({+)£  //  F2(u,  u)  du  du) 

+c,„  [»«((+)/-  //  Fi(u,  u)  du  du  +  cos((+)£  //  Fj(u,  u)  du  du] 

kLFMd'1 


(177) 


Plugging  in  expressions  for  the  derivatives  of  the  integrals,  and  implementing  Equation  129,  yields 


Z*mn  =  esncostf+)e-"*'‘ 
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h(ir  -  l)e,„  |  «*(<♦)  l «”-*W  -  e-"*» U| 
41,0  |  -«in({+) 

_*2 

4 -  c  ‘mnra 

4; 


(178) 
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Appendix  C.  Integration  in  the  Neighborhood  of  the  Singularity 


To  fill  the  impedance  matrix,  as  defined  by  Equations  172,  175,  and  178  in  Appendix  B 
we  must  evaluate  integrals  whose  integrands  have  singularities  at  the  origin.  Specifically,  these 


integrals  are 

K,(y]  =  /(_  »•■»[/■ 

»  (v  +  y  +  (koy/xi2  +  v7J  dudv 

(179) 

K,(v)  =  /£««[/ 

n  +  y  +  ^ j  e~it%uH^  {ko\J u?  +  v2j  du  dv 

(180) 

K,(y)  =  K,(y) 

(181) 

Region  Sj  is  a  circle  of  radius  R  if  the  match  point  lies  within  the  dielectric,  or  a  semicircle  of 


radius  R  if  the  match  point  lies  on  the  conductor.  Region  5j  is  illustrated  by  Figure  19.  We  are 


Figure  19.  Region  of  Integration  About  the  Singularity 
also  interested  in  evaluating  partial  derivatives  of  these  are  given  simply  by 

0 

UK, 


Bt  ""i 

±r  - 
_ 


(182) 

(183) 


04 


(134) 


We  consider  separately  the  cases  in  which  the  match  point  lies  inside  the  dielectric  or  on  the 
conductor. 


C.J  Match  Point  Inside  the  Dtclectnc 

When  the  match  point  lies  inside  the  dielectric  (y  ^  -£),  region  Si  is  a  circle  of  radius  R. 
Using  Euler’s  Identity  to  turn  the  trigonometric  terms  into  exponentials  and  switching  to  polar 
coordinates  (u  =  pcos9,v  =  p sin0)  allows  us  to  write 

_ it _ 

_ h _ 

r  r  c-i^-u'',+‘~^^ui2\t0p)pdpdo  ties) 

2 J  Jo  Jo 

+  (186) 

These  integrals  have  no  dependence  on  the  match  point  (x,  y)m  and  thus  need  only  be  evaluated 
once,  with  the  results  applied  to  all  impedance  matrix  elements  whose  match  point  lie  in  the  interior 
of  the  dielectric  We  reformulate  the  exponential  using  the  identities 


i4cos^  +  Bsin^  =  \/i42  +  B2cos^ 

>-"ctan(§)] 

(187) 

arctan  +  arctan  =  8gn(A)x 

(188) 

fA\  .  f-A\ 

arctan  f—J  +  arctan  =  0 

(189) 

/n?  +  Qn 2  =  Cr^o 

(190) 
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and  write 


/,  =  J*pHg\k0p)  ,»*/!».*■•[♦♦■■«•■(£)-»]  iedp 

h  =  /%42>(W)  ["  dp 

Jo  Jo 


(191) 

(192) 


Here  we  take  the  range  of  the  arctan  function  to  be  (-r,  *}  and  must  be  careful  to  preserve  the 
“sense”  of  the  fraction  (i  e  ,  arctan  (^f)  /  arctan  (^))-  Since  the  integration  in  0  covers  a  full 
period,  the  phase  terms  in  the  cosine  functions  don’t  matter,  and  integration  in  6  yields  a  zero-order 
Bessel  function  (17:360)  and,  rather  fortuitously,  /j  =  We  may  now  write  the  one-dimensional 
integral  in  p  as 


tR 

h  =  h  -  2x  /  M^k0p)H^\k0p)pdp 
Jo 


(193) 


We  scale  the  integration  so  that  it  is  with  respect  to  the  argument  of  the  Hankel  function  (r  =  top), 
and  write 


A  ft  (K 

/ l  =  =  Jp{p/erT)  dr 


(194) 


Requiring  y/c^koR  <  3,  or  Q  we  may  replace  the  Bessel  and  Hankel  functions  by  their 

seven-term  small- argument  polynomial  expansions  (,as  given  in  Appendix  D)  Integrating  term-by- 
term,  we  find 


j  ~  j  ..  QtCr  (*oR\  ’ 

[“‘  -  jh  (in(tofi)  -  2(75T+Tj)]  <195» 


Plugging  this  result  into  Equation  185  and  186  allows  us  to  write  the  final  result  as 


?»,  =  [/»  (v  +  !)]  A 

Kv+S]/i 


=  «* 


(196) 

(197) 
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C.2  Match  Point  on  ihc  Conductor 

If  the  match  point  (2,  y)m  lies  on  the  conductor,  then  we  define  the  region  Si  to  be  a  semicircle 
of  radius  R  centered  at  the  singular  point  (u,u)=(0,0),  extending  into  the  dielectric  In  a  manner 
identical  to  that  used  in  the  previous  section,  and  noting  that  y  +  £  =  0,  we  may  write  If  and 


_ U _ 

_ h 

1  J*  e->*/—‘+f.  “•'>/7J1>(*0pWp<» 

ftnf  =  jA  +  jA 


(198) 

(199) 


Converting  to  polar  coordinates,  we  recast  I\  and  I2  as 

A  -  j”f  HI  d<t>  dp 

A  =  £t>  ujP(tof)  dtdp 


(200) 

(201) 


We  shift  the  interval  of  integration  on  $  and  scale  the  integration  in  p  and  write 


h 

*2 


=  rj  /  [  '  ^rm,dedr 

*oJo  Arei»n({ J)-. 

*0  Jo 


(202) 

(203) 


Again  requiring  £  <  we  may  replace  the  exponential  and  Hankel  func*  ons  by  their  small- 
argument  polynomial  expansions  (aee  Appendix  D)  Integrating  term-by-term,  we  find  the  final 
expressions  for  A  and  7j  to  be 


A 


(to RY+U 
fcifcS  i!3“  '’  +  24  +  2 


X 
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R)-utd)Cl!lr'edS  (204) 


,  „  riW  WEL 

2  -  ““  i!3’*  i  +  2t  +  2  X 


{“‘-^h«)-JT^2]}/.' 


(ft)- 

uctu,(jj) 

"■“(fc)-' 


cos*  6  d$ 


(205) 


The  integral  in  8  is  straightforward  (12.132)  but  has  differing  forms  for  i  even  or  odd.  Plugging 
Equations  204  and  205  into  Equations  198  and  199,  we  may  write 


if,.,  =  Yj(h  _  h)  (206) 

If;',,  =  !<A  +  /2)  (207) 


The  exponential  series  in  Equations  204  and  205  are  taken  out  to  the  N  +  Ith  term.  The 
alternating  nature  of  these  series  allows  us  to  bound  the  error  associated  with  their  truncation  by 
the  sum  of  the  first  two  terms  neglected  (1293-4).  For  instance,  letting  JV  =  10  limits  the  relative 
truncation  error  in  Equation  204  to  less  than  10“d. 
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Appendix  D.  Some  Useful  Series  Expansions 


To  fill  the  impedance  matrix,  as  defined  by  Equation  185  and  198  in  Appendix  C,  we  must  raa- 
lytically  evaluate  integrals  that  do  not  appear  in  integral  tables  (12).  We  evaluate  these  integrals  by 
expanding  the  functions  in  the  integrand  by  their  appropriate  small-argument  finite  senes  approx¬ 
imations,  multiplying  the  series  together,  interchanging  orders  of  integration  and  summation,  and 
integrating  term-by-term  In  this  Appendix,  we  present  the  series  expansions  used  to  approximate 
the  functions  appearing  in  the  integrands  of  Equations  172,  175,  and  178  in  Appendix  B. 

D  1  Bessel  Functions  of  Order  Zero  and  One 

Abramowitz  &  Stegun  (17  369-370)  present  the  following  small-argument  polynomial  approx¬ 
imations  for  zero  and  first-order  Bessel,  Neumann,  and  Hankel  functions  The  arguments  are  real 
and  confined  to  the  interval  (0, 3). 


Jo{z) 

i=0 

(208) 

Vo(x) 

*  (;)'"(!)«■>* j>(!)’ 

(209) 

Hi%) 

=  Mz)  -  jVo(x) 

»=0 

(210) 

MCx)h£\x) 

=  £X>.C*[ot  “  !"(*))  (I)***’ 

1=0  4=0 

(211) 

The  constants  a, ,  bt%  a,,  and  /?,  are  given  in  Table  4  The  absolute  errors  associated  with  these 
approximations  are  less  than  5  x  10~8. 

Similarly,  the  small-argument  polynomial  expansions  for  Bessel,  Neumann,  and  Hankel  func¬ 
tions  of  order  one  are 

s  2i 

Hz)  -  *£«■(§)  (212) 

{=0 
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i 

a. 

i>. 

<*i 

0. 

0 

1.0000000 

0  36746691 

1. 0000000+ j 0.0738043 

0  63661977 

i 

-2.2499997 

0.60559366 

-2. 2499997- j 1.6984537 

-1.43239430 

2 

1.2656208 

-0.74350384 

1. 2656208+ j 1.3019858 

0  80571923 

3 

-0.3163866 

0  25300117 

0. 3163866- j 0.3926136 

-0.20141797 

4 

0.0444479 

-0.04261214 

0 . 0444479+ jO . 0622267 

0  02829641 

5 

-0.0039444 

0.C0427916 

-0. 0039444- jO. 0060197 

-0.00251108 

6 

0  0002100 

-0.00024846 

0.0002100+j0. 0003411 

0.00013369 

Table  5.  Coefficients  for  Bessel  and  Neumann  Function  Expansions,  Order  One 


i 

Cx 

d, 

Tl 

0.50000000 

-0  6366198 

l 

-0.56249985 

0.2212091 

2 

0.21093573 

2.1682709 

3 

-0.03954289 

-1.3164827 

4 

0.00443319 

0  3123951 

5 

-0.00031761 

-0  0400976 

6 

0.00001109 

0.0027873 

tt(z) 

»?\*) 


JiW-yv.w 


(213) 

(2M) 


where  the  constants  c»  and  d,  are  given  in  Table  5  The  absolute  errors  associated  with  these 
approximations  are  less  than  1  1  x  10-7. 


D  2  Hankel  Functions  MuUtpUcd  by  Ezponenitals 

We  need  to  evaluate  integrals  of  Hankel  functions  multiplied  by  complex  exponentials,  r  task 
most  easily  performed  by  small-argument  polynomial  approximations  Hankel  function  polynomial 
approximations  are  given  in  the  previous  section,  exponentials  may  be  approximated  by  (12  22,34) 


<?Cx 


1*1  <°° 


(215) 
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c<*(Cx)  ~  (216) 

The  exponential  series  is  written  in  this  form  to  highlight  its  alternating  nature,  which  allows  us 
to  bound  the  error  associated  with  truncating  the  infinite  series  by  the  magnitude  of  the  first  term 
neglected  (1:293-4).  Combining  series,  we  can  then  write 

[> + (217> 

1=0  i  =  0  '  '  L  J 

cos (Cx)H<?\z)  =  (218) 

1=0 1=0  *  ' 
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Appendix  E.  Derivatives  of  the  Two-Dimensional  Green’s  Function 


The  two-dimensional  lree-space  Green’s  function  is  given  by  (3  902) 


G(x,  »!*',✓)  =  jjHpihR) 

R  =  \/(t  -  z')2  +  (v  -  l/y 


(219) 

(220) 


where  f/J2'  is  the  Hankel  function  of  the  second  kind  of  order  0.  We  wish  to  find  expressions  for  t 
ef  ’  an<*  =  start  with  e?  Using  the  chain  rule  and  a  Hankel  function 

recurrence  relation  (17  361),  (koR)  =  —H^\k0 R),  we  find  the  first  partial  derivative  with 
respect  to  x  to  be 


dG  _  dGdR 
dz  ~  dR8z 


=  j-°-^H?\k0R) 


(221) 


where  the  prime  on  the  Hankel  function  denotes  differentiation  with  respect  to  argument  By 
symmetry, 


8G  ,io(sr  —  lO  m 
r~h'  {t°R) 


The  second  derivative  with  respect  to  x  is 


(222) 


d7G  _  d_8G 
dz2  ~  dr  dx 

Using  the  product  rule,  chain  rule,  and  another  Hankel  function  recurrence  relation  (17:361), 
4%oR)  -  j^//[J)(hoR)  =  -H?\koR) 
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we  find 


d7G 

3x2 


=  f  +  *>S  35  (5) 


(x  —  xf)  5/2  d 


/?  5x  d/2 


■M'HttR)} 


=  f  |^£«)  _  ^*<2>(ioK)  +  ^2£flp/(WD} 


J*Q  I  gl!Wg>  _  (*-«')*  „m 


T\— 5-  ~  Tw'  (ioR) 


Of  special  interest  is  the  term 


5JG 
3x 2 


[*o  +  j  G  =  *oG  + 

=  |«J5)(iofl)  -  f  -  — /jr^2,(*oK)} 


=  ~  {4v(toR)  -  ^h[7\kr)  + 


By  symmetry, 


0  = 


[^01c  =  UiH<:  ihn>  - -■^■^){toR)} 


Examination  of  these  equations  shows  that 


V’G  +  i§G  =  j-f  +  0  +  ij  j  =  0 


(223) 


(224) 


(225) 

(226) 


(227) 
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Using  the  recurrence  relation  h\2*  {koR)  -  —  -H^ikoR)  (17  361),  we  find 

£&  to  be  given  by 


jko(t-z')dR  d  { Hf'ikcR)] 

4  By  dR  \  R  ) 

jkc(z  -  z')(y  -  ]/)  k0 
-  - -  R  \rfl  ^ 

k%  (z  -  x')(y  - !/) 


(koR)  - 


H^lhR)) 
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